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1 .0  INTRODUCTION 


This  report  summarizes  work  done  in  connection  with  the  high- 
resolution  transmission  molecular  absorption  database  HITRAN^.  This 
database  is  a  compilation  of  spectroscopic  parameters  from  which  one  is 
able  to  calculate  and  predict  the  transmission  and  emission  of  radiation  in 
the  atmosphere.  The  work  presented  here  spans  a  wide  range  of  the 
problems  one  encounters  when  maintaining  a  database.  In  this  report  we 
discuss  calculations  of  fundamental  parameters  needed  by  the  database, 
calculations  of  data  needed  for  application  of  the  database  to  specific 
problems,  the  creation  of  maintenance  and  updating  programs  which  help 
eliminate  the  many  errors  that  can  occur  when  dealing  with  large  amounts 
of  data,  creating  the  newest  version  of  the  database,  and  we  discuss  the 
overall  philosophy  used  in  the  development  of  the  1991  database. 

In  the  second  section,  the  work  done  in  computing  total  internal 
partition  sums  for  temperatures  from  70  to  3000  K  for  the  species  on  the 
HITRAN  database  and  the  fitting  of  this  data  to  particular  models  is 
discussed.  The  total  internal  partition  sums  (TIPS)  allows  the  application 
of  the  database  to  other  temperatures  as  well  as  the  calculation  of  the 
transition  moment  squared  or  Einstein  coefficients.  Section  3  gives  the 
programming  philosophy  adopted  for  all  programs  that  interact  with  the 
database.  Various  programs  used  for  updating,  maintenance,  and 
producing  the  new  versions  are  presented  and  the  interrelationships 
between  the  programs  shown.  The  fourth  section  describes  the  effort 
behind  the  1991  version  of  the  HITRAN  database  describing  all  differences 
one  will  observe  between  the  1986  version  and  the  1991  version  as  well 
as  the  structural  changes  made  to  the  database.  Section  5  discus.ses  the 
calculation  of  pressure-broadening  coefficients  for  various  species  on  the 
database.  The  application  of  direct  access  files  to  create  mixed  databases 
of  experimental  and  theoretical  pressure-broadening  coefficients  for 
addition  to  HITRAN  is  discussed.  Section  6  gives  the  relationships  between 
the  line  intensities,  transition  squared,  and  the  Einstein  coefficient.  Section 
7  reviews  the  two  workshops  that  were  held  on  the  HITRAN  database  and 
gives  some  of  the  views  of  the  author.  Finally  section  8  gives  some  of  the 
presentations  and  publications  derived  from  this  work. 
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2.0 


TOTAL  INTERNAL  PARTITION  SUMS 


During  the  period  of  this  contract  the  total  internal  partition  sums 
(TIPS)  were  calculated  for  most  of  the  isotopic  variants  on  the  database. 
This  augmented  some  earlier  work  done  for  the  principal  atmospheric 
species.  The  total  internal  partition  sum  is  given  by 

Q(T)  =  X  (2J+1)  (2S+1)  (2-8x0)  (gn)  e-E/kT 

all  levels  ,  ( 1 ) 

where 

(2J+1)  J  spatial  degeneracy 

(2S+1)  ms  spatial  degeneracy 

(2-8xo)  lamda  doubling  degeneracy 

(gn)  nuclear  spin  (state  dependent  and  state  independent) 

degeneracy 

E  Total  energy  of  the  level,  electronic,  vibrational, 
and  rotational. 

The  original  intent  was  to  extract  quantum  numbers  and  energies  directly 
from  the  database  and  perform  the  sum  in  equation  (1).  It  was  realized 
that  for  many  of  the  species  on  the  database  the  resulting  TIPS  would  not 
be  converged  at  the  elevated  and  even  intermediate  temperatures  under 
consideration.  Thus  it  was  decided  to  calculate  near  complete  sets  of 
energy  levels  via  the  Hamiltonian  of  the  species  in  question.  Attempts 
would  be  made  to  go  to  high  enough  energy  such  that  these  terms  in  the 
sum  lEq.  (1)]  would  have  little  to  no  effect  on  the  final  results. 

Direct  summation  of  equation  (1)  over  all  vibration-rotation  levels 
was  done  whenever  possible.  For  many  of  the  species  however,  molecular 
constants  are  not  available  for  all  vibrational  bands.  The  TIPS  via  a  direct 
sum  for  this  situation  would  have  many  energy  levels  missing  and/or 
would  have  a  cutoff  in  energy  that  is  too  low.  For  these  cases  it  was 
assumed  that  the  vibrational  and  rotational  energies  are  separable  and  the 
partition  sum  can  be  written  as  a  product  of  the  vibrational  partition  sum 
and  the  ground  state  rotational  partition  sum, 

Q(T)  =  Qvib(T)xQ,odT) ,  (2) 
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The  effect  of  a  cutoff,  Figure  la,  is  to  remove  all  energies  above  the  cutoff 
value  from  entering  into  the  sum.  Applying  a  cutoff  to  the  ground  state 
and  then  applying  the  product  approximation,  Qvib(T)*Qrot(T),  has  the 
effect  of  including  many  energy  levels  above  the  ground  state  cutoff  as 
shown  in  Figure  lb. 


Figure  la.  Energy  cutoff  applied  to  ro-vibrational  energy  levels. 


Vibrational  Levels 


Figure  lb.  Energy  cutoff  applied  to  ground  state  energy  levels, 
followed  by  the  product  approximation. 


When  the  product  approximation  was  made  both  Qvib  and  Qrot  were 
calculated  by  direct  summation.  In  the  calculations,  the  energies  of  the 
ground  vibrational  states  were  determined  from  the  Hamiltonian  for  each 
species  in  question.  The  vibrational  partition  sums  were  calculated  using 
the  vibrational  energies  in  the  literature  or  in  some  cases  estimated  for 
lesser  isotopic  species.  The  vibrational  sum  is  usually  over  5  to  10  terms 
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and  the  rotational  sum  can  be  over  several  hundred  to  several  thousand 
levels.  The  product  approximation  was  shown-  to  be  accurate  to  a  few 
tenths  of  a  percent  for  water  vapor,  carbon  dioxide  and  ozone.  When 
energy  levels  were  not  available  for  direct  summation  then  the  classical 
values  were  used. 

The  TIPS  were  to  be  calculated  at  a  number  of  temperatures  and 
these  fit  by  different  expressions.  After  some  consideration  we  settled  on 
fitting  in  three  temperature  ranges;  70-415K  in  5K  steps,  385-2005K  in 
20K  steps,  and  by  extrapolation  of  the  1900  and  2005K  points  to  the  2005- 
3000K  rang*.. 

Several  representations  for  the  temperature  dependence  of  the  total 
internal  partition  sums  were  tried  by  fitting  to  the  data.  The  first 
expression  used  came  from  the  work  of  McDowell^  and  had  the  form 

where  P  =  hcB/kT  (B  =  rotational  constant)  and  o*  =  (21y  +  l)n/a  (ly  = 
nuclear  spin  of  the  Y  nuclei,  o  =  classical  symmetry  number)  and  D  = 
centrifugal  distortion  constant.  In  this  expression  a,  p,  D  and  B  were 
treated  as  adjustable  constants. 

A  simple  polynomial  expression  in  temperature  with  four  adjustable 
constants, 

2  3 

Q(T)=a+bT+cT  +dT  ^  (4) 

was  then  studied. 

There  is  another  form  that  is  sometimes  used.  This  considers  the 
natural  logarithm  of  Q(T)  vs,  1 /temperature.  Several  points  are  stored  and 
values  at  other  temperature  can  be  obtained  from  the  formula 

ln(Q(T))=a,b 

The  constants  a  and  b  are  obtained  from  the  stored  points.  For  example  S. 
Clough  uses  this  method  in  the  FASCODE^  program  and  H.  Pickett  uses  it  in 
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the  JPL  Catalogue^. 

Test  were  performed  by  fitting  expressions  (3)  through  (5)  to  Q(T) 
data  and  then  recalculating  Q(T)  and  comparing  this  to  the  original  QC'D 
data.  The  results  show  that  the  polynomial  form,  Eq.  (4),  for  the  partition 
sum  gave  the  most  accurate  results  and  is  the  simplest  for  recalculation  of 
Q(T).  This  was  the  form  that  was  chosen  to  represent  the  partition  sums. 

The  constants  in  equation  (4),  the  a,  b,  c,  and  d  coefficients,  were 
obtained  by  least-squares  min-max  fitting  procedure.  The  fit  in  the  70- 
41 5K  range  was  constrained  to  better  than  1%  error.  The  fitting  errors  are 
always  better  than  1/2%.  The  385-2005K  range  was  constrained  to  1% 
error  at  4I5K  and  allowed  to  float  to  5%  at  2005K.  This  is  because  the 
partition  functions  are  less  accurate  at  the  high  temperatures.  This  fitting 
procedure  dO'  -•  not  jeopardize  the  fit  at  low  temperatures  where  Q(T)  is 
accurate.  The  final  high  temperature  range  results  vary  from  1/2%  at 
41 5K  to  4%  at  2005K  as  the  worst  case. 

For  the  high  temperature  extrapolation  of  the  partition  functions  in 
the  range  2005-3000K,  the  linear  parts  of  several  forms  were  considered. 
The  linear  part  was  chosen  because  of  the  poor  extrapolation  qualities  of 
polynomial  representations.  The  forms  chosen  were  the  form  used  by  S. 
Clough'*  and  H.  Pickett^ 

lnQ(T)  =  ^.a„ 

from  Irwin^ 

In  Q(T)  =  ai  ln(T)+ ao  (7) 

and  two  classical  (rotational)  forms 

Q(T)  =  aiT  +  ao  and  In  Q(T)  =  ai  T  + ao .  (8) 

To  test  the  methods,  the  total  internal  partition  sums  were  calculated 
for  H35C1  from  T  =  2000  K  to  3000  K  in  100  K  steps.  For  this  species  the 
energies  were  calculated  to  44,500  cm- 1  giving  the  TIPS  an  uncertainty  of 
less  than  2%.  Using  the  coefficients  previously  obtained,  the  extrapolated 
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partition  sums  were  calculated  for  each  method  and  the  %  error  with  the 
direct  sum  computed.  These  values  were  plotted  in  Figure  2.  Similar 
results  were  obtained  for  CO2  ^nd  CO.  All  results  indicate  that  expression 
(2)  gives  the  best  extrapolation  of  the  methods  tested. 


The  remaining  terms  needed  to  evaluate  equation  (1)  (once  the 
energies  are  obtained)  are  the  statistical  fr'ctors.  Some  of  these  are 
dictated  by  the  quantum  numbers,  thus  they  are  state  dependent.  Others 
are  independent  of  the  particular  vibrational-rotational  state  and  are  state 
independent  factors  and  can  be  brought  out  of  the  summation  in  Eq.  (1). 
Because  these  state  independent  factors  are  constant  for  a  given  species 
and  usually  ratios  of  the  partition  functions  are  used  they  arc  commonly 
omitted.  The  internal  partition  sums  calculated  in  this  work  contain  all 
terms  explicitly.  Thus  they  can  be  used  to  calculate  thermodynamic 
functions,  etc. 

This  includes  nuclear  statistical  factors  that  are  state  independent. 
These  terms  arise  from  the  nuclear  statistics  of  nuclei  that  do  not  couple 
with  the  rotational  wavefunctions  because  of  symmetry.  The  nuclear 
statistical  factor  generates  additional  multiplying  terms  of  gj  =  (2I-I-1)  for 

each  nuclei  that  is  not  coupled  with  the  rotational  wavefunctions,  where  1 
is  the  spin  of  the  nucleus  in  question.  These  have  been  calculated  and  are 
given  in  Table  1. 
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Table  1.  State  Independent  Nu:lear  Spin  Degeneracy  Factors 


AFGL  AFGL 

isetope  gj  =  n(2Ij  +1)  Isotope 

Molecule  Code  Molecule  Code 


H20 

161 

1 

HNO3 

146 

181 

1 

171 

6 

OH 

61 

162 

6 

81 

62 

C02 

626 

1 

636 

2 

HF 

19 

628 

1 

627 

6 

HQ 

15 

638 

2 

17 

637 

12 

828 

1 

HBr 

19 

728 

6 

1 1 

03 

666 

1 

HI 

17 

668 

1 

686 

1 

ao 

56 

76 

N20 

446 

9 

456 

6 

OCS 

622 

546 

6 

624 

448 

9 

632 

447 

54 

822 

CD 

26 

1 

H2OO 

26 

36 

2 

136 

28 

1 

128 

27 

6 

Hoa 

165 

CH4 

211 

1 

167 

311 

2 

212 

3 

N2 

44 

O2 

66 

1  * 

HCN 

124 

68 

1 

134 

67 

6 

125 

gj  =  n(2Ij+l) 


6 

2 

2 
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8 
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12 
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1 
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1 

1 

2 

1 

8 

8 
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Tabic  1 

continued 

fO 

46 

3 

56 

2 

48 

3 

Sp2 

626 

1 

646 

1 

NO2 

646 

3 

NH3 

4111 

3 

5111 

2 

0130 

215 

4 

217 

4 

H2O2 

1661 

4 

C2H2 

1221 

1* 

1231 

8 

C2H6 

1221 

64 

PH3 

nil 

2 

*  not  applicable  due  to  symmetry 


2.1  Water  Vapor,  H2O 

For  water  vapor  it  is  not  possible  to  calculate  all  vibration-rotation 
ene»'gies  beyond  =  7000  cm'^,  so  we  have  calculated  the  total  internal 
partition  sums  using  the  approximation  Q(T)  =  Qvib(T)*Qrot(T).  This  has 
the  effect  of  approximately  including  many  of  the  high  lying  energy  states 
that  are  removed  when  a  cutoff  is  applied.  The  rotational  partition  sums 
were  calculated  with  the  Watson  Hamiltonian  constants  of  Flaud  et  al.^  for 
H2*^0,  H2^^0,  and  H2^^0.  The  vibrational  energies  were  taken  from 
HITRAN  to  do  a  direct  sum  for  the  vibrational  partition  sum. 

2.2  Carbon  Dioxide,  CO2 

For  *  20160 2  the  molecular  constants  of  Wattson  and  Rothman^  from 
the  direct  numerical  diagonalization  (DND)  method  were  taken  to  generate 
the  vibrational-rotational  energy  levels.  This  produced  a  set  of  energies 
complete  to  16000  cm-l  with  a  cutoff  of  20000  cm-l.  For  1301602  and 
160120180,  the  constants  of  Rothman  et  al.9  were  used,  this  generated 
energies  up  to  14055  cm'l  and  13860  cm-l  respectively.  For  these 
species  direct  summations  were  used  to  compute  the  TIPS. 

For  the  other  five  less  abundant  species  the  energy  levels  were  not 
available  for  all  vibrational  bands.  For  these  species,  1^0120170,  16ol30l8o, 
18ol20l8o,  17ol20l8o,  the  product  approximation  was  employed 
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to  evaluate  the  partition  sums.  The  ground  rotational  parameters  are  from 
Chedin  and  Teffo^O,  The  vibrational  partition  sums  used  the  vibrational 
energies  reported  in  Reference  1.  The  product  Qvib  *  Qrot  was  done  to 
give  the  total  internal  partition  sums  for  these  species  of  C02- 

2.3  Carbon  Monoxide,  CO 

The  energy  levels  of  13(^160,  12(^180,  l2(;;17o,  and  *3cl8o  were 

calculated  using  the  Dunham  expansion  with  the  constants  of  Guelachvili 
et.  al.ll.  The  energies  were  calculated  for  v  =  0  to  9  and  J  =  0  to  74.  The 
range  of  v  and  J  considered  generated  a  complete  set  of  energy  levels  to 
10000.  cm"^  and  we  expect  the  total  partition  sums  will  converge  for  all 
temperatures. 

2.4  Monodeuteraied  Methane 

For  CH3D  the  energy  levels  have  not  been  calculated  yet.  To  compute 
the  total  internal  partition  sum  for  this  species  we  have  taken  the 
rotational  partition  sum  from  Robiette  and  Dang-Nhu'^  at  300  K, 
Qrot=9095.0,  and  computed  the  temperature  variation  as 

(Q(T))  =  Q(300K)x(;^P 

'300'  .  (9) 

The  vibrational  partition  sum  was  evaluated  by  direct  summation  of  the 
vibrational  energies  given  in  HlTRANl.  The  total  internal  partition  sum 
was  then  computed  via  Qvib'Qrot- 

2.5  Ammonia,  NH3 

The  energy  expression  for  ammonia  is  given  reference  13  with  the 
inversion  energy  was  calculated  using  the  Fade  approximation  of  Young 
and  Young '4  The  constants  used  for  14nh3  are  from  Poynter  and 
Margolisl3  with  corrections  to  the  constants  C,  Dk,  Hk,  and  Lk  from  Cohen 
et.al.l5  The  inversion  energy  constants  were  taken  from  reference  14. 
Symmetry  rules  for  the  allowed  energy  levels  from  Dowlingl^  were  used  to 
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write  an  energy  level  program.  With  the  program,  energies  were 

calculated  for  J  -  0  to  50  (=22000  cm'^ 

The  statistical  factor,  g^,  for  explicit  rotational-inversion  levels  was 
taken  from  Bunkerl^.  With  these  factors  the  rotational  partition  sums 
were  calculated  for  the  temperatures  in  question. 

For  the  symmetric  and  antisymmetric  rotation-inversion 

constants  were  taken  from  Carlotti  et.  al.*^  with  vq  for  inversion  from 
Sasadal^.  Energies  were  calculated  for  J  =  0  to  50  (=22000  cm-l). 

Vibrational  energies  were  taken  from  subroutine  MOLEC^O  (part  of 
FASCODE^)  to  give  the  total  internal  partition  sums  for  these  species. 

2.6  Sulfur  Dioxide,  SO2 

The  energies  of  the  ground  vibrational  state  were  calculated  using 
the  Watson  Hamiltonian  constants  of  Carlotti  et  al.^l  for  the  species  ^25100 2 
and  the  constants  of  Pine  et  al.^^  for  the  species  ^"^8*^02.  Calculations  were 
done  for  J  =  0  to  95  and  levels  up  to  =7000  cm'^  were  retained.  The 
vibrational  partition  sum  used  the  vibrational  energies  from  Reference  22. 
These  are  for  the  principal  isotopic  species  only  thus  the  sums  for  the 
3451602  species  will  be  more  in  error,  however,  this  is  expected  to  be  very 
small.  The  rotational  and  vibrational  partition  sums  were  calculated  and 
the  product  taken  for  the  TIPS 

2.7  Nitric  Acid,  HNO3 

For  nitric  acid,  the  Watson  Hamiltonian  constants  that  were  used 
came  from  the  work  of  Maki  and  Wells23.  The  calculation  of  the  rotational 
energies  considered  J  =  0  to  103  to  produce  some  10816  energy  levels  to  = 
4500  enr  1 .  The  vibrational  partition  sum  employed  the  vibrational 
energies  published  in  the  HITRAN  paper*  as  compiled  by  Rinsland-^.  The 
total  internal  partition  sum  was  calculated  by  the  approximation, 
Qvib*Qrot. 
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2.8  Hydroxyl  Radical,  OH 


The  calculation  of  the  energy  levels  followed  the  methodology  of 
reference  25  for  ^^qh,  and  reference  26  for  and  l^OH.  The  methods 

are  a  direct  diagonal ization  of  the  three  interacting  states.  The  energies 
obtained  from  diagonalization  were  assigned  as  follows,  the  lowest  belongs 
to  the  2n3/2  state,  the  slightly  higher  value  belongs  to  the  ^ni/2  state,  and 
the  largest  value  is  for  the  state. 

A  program  was  written  to  determine  the  energies  by  diagonalizing 
the  matrices  as  a  function  of  J,  and  F  values  for  via  the  formalism 

above^S.  The  molecular  constants  used  were  from  Goldman27  for  the  Hrot 
and  Hso  terms,  from  Coxon^S  for  the  hyperfine  terms,  and  from 
Destombes25  for  the  second  order  hyperfine  terms.  Energies  E(J,  F)were 
calculated  for  all  F  states  of  a  given  J  for  J  =  1/2  to  74.5.  This  generated 
energies  up  to  6  xlO^. 

For  the  other  species,  '^qh  and  ^^OD,  constants  were  not  available 
from  the  literature  for  the  same  analysis.  For  these  species  the  formalism 
of  Beaudet  and  Poynter^^  was  used.  The  formalism  is  sirr.ilar  to  that  of 
Destombes  et.  al.  except  all  second  order  terms  and  higher  are  omitted. 
The  Calculations  were  done  for  I^OD  and  l^OH  for  J  equal  to  1/2  to  74.5  (all 
F  states).  The  rotational  partition  sums  were  then  calculated  for  each 
species.  The  vibrational  energies  were  calculated  using 

Gv  =  Gi^  + We(v'+i)-W^e(v'+^]^ 

Constants  for  l^OH  are  from  Coxon29  and  for  from  Coxon^O.  The 

constants  for  were  assumed  the  same  as  for  ^^OH.  The  total  internal 

partition  sums  were  then  calculated  via  the  approximation  Qtot  =  Qrot  * 
Qvib- 

2.9  Hydrogen  Bromide,  HBr 

The  energy  levels  of  this  diatomic  molecule  were  calculated  using  the 
Dunham  expansion  with  the  constants  of  Stocker  and  Goldman^l.  The 
energies  were  calculated  for  v  =  0  to  9  and  J  =  0  to  74  which  generates 


energy  levels  to  43811.  cni"^  thus  the  total  partition  sums  will  converge 
for  all  temperatures. 

2.10  Hydrogen  Iodide,  HI 

For  HI,  the  energies  were  calculated  by  the  Dunham  expression  using 
constants  provided  by  R.  Tipping32.  The  energies  were  calculated  for  v  = 
0  to  9  and  J  =  0  to  74  which  generates  energy  levels  to  31000.  cm*^  and 
the  total  partition  sums  will  converge  for  all  temperatures. 

2.11  Chlorine  Monoxide,  CIO 

Chlorine  monoxide  is  an  inverted  spin  doublet  with  the  ^113/2  state 
being  lower  in  energy  than  the  ^  FI  1/2  state.  the  energy  levels  were 
calculated  by  an  effective  Hamiltonian  following  Cohen  et.  al.33.  The 
molecular  constants  for  the  rotational,  hyperfine,  and  quadrupole 
interactions  were  taken  from  references  33-36.  The  rotational  partition 
sums  were  calculated  from  these  energy  levels.  The  vibrational  energies 
from  reference  36  were  taken  to  compute  Qvib.  The  total  internal 
partition  sums  were  then  calculated  by  the  product  approximation. 

2.12  Carbonyl  Sulfide,  OCS 

For  OCS  the  rotational  Hamiltonian  constants  are  from  Lovas37. 
Energy  levels  were  calculated  for  each  species  from  J  =  0  to  200,  this  yields 
energies  to  «  8000  cm"^.  The  rotational  partition  sums  were  calculated 
from  these  energy  levels.  The  vibrational  partition  sums  were  calculated 
using  the  vibrational  energies  from  Ref  1.  The  total  partition  sums  were 
done  via  the  approximation  Qtot  =  Qvib  *  Qrot- 

2.13  Formaldehyde,  H2CO 

The  rotational  energies  of  the  two  species  and  H2'3Ci<^0  were 

taken  from  Winnewisser  et  al.38  Calculations  considered  levels  from  J  =  0 
to  75  and  a  7000  cm'^  cutoff  was  applied  to  give  energies  for  the 


rotational  partition  sums.  The  rotational  partition  sums  were  then  done 
for  these  two  species.  Watson  constants  have  not  yet  been  obtained  for 
the  H2’2C*80  species.  The  vibrational  energies  for  the  principal  species  are 
from  reference  1.  These  were  used  to  compute  the  total  partition  sums  for 
this  species. 

The  total  internal  partition  sum  for  the  species  was 

calculated  via  the  assumption 

.  (II) 

The  reason  for  this  assumption  was  that  the  vibrational  energies  of  the 
Hi 3c  160  species  were  not  known.  For  H2^2C180  there  are  no  Watson 
Hamiltonian  constants  available  for  this  isotopic  species.  For  the  time 
being  the  partition  sum  for  this  species  has  been  assumed  to  equal  that  of 
the  principal  species  H2^^C160. 

2.14  Hydrogen  Hypochlorite,  HOCI 

For  the  vibrational  energies,  the  HiTRAN  values  were  used  for  both 
isotopic  species.  The  energy  levels  of  HOCI  were  calculated  using  the 
Watson  Hamiltonian  constants  of  Wells,  Sams,  and  Lafferty^^.  Some  5776 
energy  levels  were  generated  for  each  isotopic  species  (HO^^CI  and  HO^^Cl) 
by  calculating  levels  from  J  =  0  to  75  and  cutting  off  at  7000  cm"^  to  give 
2326  H03-'’C1  lines  and  2332  HO^^CI  lines.  The  rotational  partition  sums 
were  done  by  summing  over  these  energj  values  with  appropriate 
degeneracies.  The  final  total  sums  were  then  calculated  by  the  product 
formula. 

2.15  Nitrogen,  N2 

For  the  nitrogen  molecule,  the  energies  were  calculated  for  all  states 
from  v  =  0  to  9  and  J  =  0  to  75  via  the  Dunham  expansion  using  the 
constants  of  Reuter  et.  al.^9  xfijs  g^yc  all  energy  levels  from  0  to  -31000 
cm-i.  The  nuclear  statistical  factors  came  from  Hertzberg^l.  For  •"lN2  this 
gives  gs  -  6  for  i.he  symmetric  levels  (J  even)  and  gs  =  3  for  the  anti- 


symmetric  levels  (J  odd).  The  total  internal  partition  sum  was  calculated 
by  direct  summation  over  all  energies., 

2.16  Hydrogen  Cyanide,  HCN 

The  constants  for  the  three  isotopic  species  considered  for  this 
molecule  are  from  the  work  of  Maki42  with  some  updates  from  Lovas^2 
Calculation  of  levels  from  J  =  0  to  75  were  done  for  each  species.  This 
produced  energies  up  to  =  8000  cm"^.  The  rotational  partition  sums  were 
calculated  at  the  150  temperatures  using  these  energies.  The  vibrational 
energy  levels  are  available  for  from  HITRAN,  for  from 

Alpert,  Mantz,  and  Rao'^^,  and  for  from  reference  43  and  scaling  of 

the  vibrational  levels.  The  final  TIPS  were  calculated  via  the 

product  of  Qvib*Qrot- 

2.17  Hydrogen  Peroxide,  H2O2 

Hydrogen  peroxide  is  a  light  asymmetric  rotor  and  it  has  substantial 
centrifugal  distortion  contributions  to  its  rotational-torsional  spectrum. 
The  torsional  states  are  fourfold  degenerate,  in  the  lowest  torsional  state 

(n=0)  the  fourfold  degeneracy  is  broken  and  a  doublet  of  doublets  occurs. 
The  four  sublevels  are  designated  by  x  =  1,  2,  3,  4  in  order  of  increasing 
energy.  The  x  =  1  and  x  =  2  states  and  the  x  =  3  and  x  =  4  states  are 
degenerate.  The  energy  levels  were  calculated  via  ASMROT  using  the 
rotational-torsional  constants  of  Helminger  et.  al."^^  for  the  x=l,2  and  x=3,4 
torsional  states.  Energies  were  calculated  for  J=0  to  35  (=  10000  cnr*)  and 
were  used  to  calculate  the  ground  state  partition  sum.  The  vibrational 

energies  from  MOLEC  were  used  to  produce  the  vibrational  partition  sums 

and  the  product  taken  for  the  total  internal  partition  sums. 

2.18  Acetylene,  C2H2 

For  acetylene  the  molecular  parameters  for  the  two  species  are  from 
Hietanen  et.  al.'^-*’’^^’.  Energies  were  calculated  up  to  J  =  100  (12,000.  cm-1) 
for  the  ground  state  of  each  species  The  ground  state  partition  sums  were 


evalu?>ted  and  for  the  principal  species  agree  very  well  with  the  value 
given  by  McDowell^^  (at  100  K,  Q(here)  =  118.8487  and  Q(McDowell)  = 
118.8263).  The  vibrational  energies  were  from  MOLEC.  The  total  internal 
partition  sums  were  done  at  the  150  temperatures  of  the  grid  used  for 
fitting. 

2.19  Phosphine,  PH3 

The  energy  expression  for  PH3  is  given  reference  48  .The  inversion 
energy  is  negligible  for  PH3  and  was  omitted.  The  molecular  constants  for 
PH3  were  taken  from  Maki  et.  al.^^.  Energies  were  calculated  for  J  =0  to  75 
(«25000  cm-*)  from  which  the  ground  state  partition  sums  were  evaluated 
at  the  desired  temperatures.  The  statistical  factors,  gj's,  are  given  in 
reference  49.  Vibrational  energies  were  taken  from  reference  1  to  give 
the  total  internal  partition  sums. 


3.0  PROGRAMMING 

There  were  many  programming  tasks  during  this  contract  period. 
Many  of  the  programs  created  were  for  one  time  use  and  were  stand  alone 
programs.  However,  there  are  many  programs  that  are  used  often  and 
interact  with  the  evolution  of  the  database.  Some  of  these  are 
maintenance  programs  for  the  database  and  others  perform  operation  on 
the  database  and/or  auxiliary  files.  In  order  to  create  a  set  of  maintenance 
programs  that  will  remain  compatible  with  the  database  now  and  in  the 
future,  a  modular  programming  structure  was  adopted.  The  reason  for  this 
is  that  many  of  the  programs  use  subroutines  and  data  that  are  common  to 
all,  e.g.  the  vibrational  state  indices.  If  this  information  is  updated  in  one 
program  it  must  then  be  updated  in  all  programs  calling  on  that  data. 
Unfortunately  the  editing  of  all  programs  containing  the  data  is  often 
overlooked.  With  the  modular  structure  common  information  is  kept  in 
individual  files  consisting  of  one  to  several  subroutines  and  a  BLOCK  DATA 
subroutine.  When  changes  are  made  to  the  module,  the  module  is 
compiled  and  the  object  code  linked  with  the  appropriate  maintenance  or 
operational  programs.  Figure  3  shows  some  of  the  interrelationships 


between  the  modules  and  the  maintenance  and  updating  programs. 
(Unfortunately,  other  programs  may  also  possess  a  dependence  on  these 
modules.  An  example  is  GENERATOR.FOR  which  creates  the  final  set  of  CO2 
parameters  and  must  be  consistent  with  the  vibrational  indices  chosen  for 
the  current  edition  of  HITRAN.) 
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Figure  3.  Interrelationships  between  modules  and  programs. 


3.1  Program  SELECT 

SELECT  is  the  u.ser-friendly  interface  to  the  HITRAN  database.  Its 
function  is  to  allow  users  to  remove  from  the  database  the  portion  of  the 
transitions  they  require.  Many  changes  have  been  made  to  the  program 
SELECT  to  produce  the  new  version  SELECT91.  These  include;  the  coupling 
of  a  wavenumber/block  correspondence  table  to  decrease  the  time 
necessary  to  step  ahead  to  the  beginning  wavenumber  of  the  requested 
search;  temperature  dependent  output  files  for  a  range  of  temperatures 
from  75  to  2005  K;  a  faster  .search  procedure;  an  option  to  create  batch  run 
input  file  in  the  user-friendly  mode, and  input  to  the  program  has  also 
been  made  easier  and  more  fool-proof. 

The  wavenuniber-block  correspondence  table  has  been  used  to 
improve  speed  of  selection  to  up  to  a  factor  of  -9  compared  with  the  1982 
version.  This  is  accomplished  by  the  new  file  structure  shown  in  figure  4. 
This  structure  allows  very  fast  selection  for  all  molecule,  particular 


molecule/all  isotopes,  and  particular  molecule/particular  isotope  searches. 
The  increased  speed  is  accomplished  by  writing  a  1986  format  file  of  the 
selected  transitions  and  using  the  wavenumber-block  correspondence 
table  to  quickly  remove  the  transitions  and  encoding/decoding  the 

niinimum  number  of  bytes  in  the  data.  This  is  shown  in  figure  5. 

The  time  consuming  encoding/decoding  is  only  done  when  a  user 
request  other  than  the  1986  standard  format.  These  formats  are  shown  in 

figure  4  and  constitute  what  are  labeled  file  3  (hard  copy  output)  and  file 

4  (line  file  in  another  format  or  ai  a  temperature  other  than  296  K).  The 
encoding/decoding  is  only  done  on  the  selected  lines  saving  a  great  deal  of 
time. 

The  1991  version  of  HITRAN  has  the  wavenumber-block  table  in  the 
first  two  blocks  along  with  a  version  number  that  will  link  particular 
versions  of  HITRAN  and  SELECT.  All  changes  will  be  downward 

compatible.  This  also  implies  that  mi.smatching  of  SELECT  and  HITRAN 
versions  will  default  to  the  slow  selection  mode. 
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File  structure  for  SELECT  1991 
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Figure  5. 


Fast  read  structure  employing  the  wavenumber-block 
correspondence  table. 


The  temperature  dependent  output  option  allows  a  user  to  generate 
a  file  of  spectral  data  at  temperatures  other  than  296°  K.  The  current 
limitations  are  T  =  70  to  3000°  K.  The  corrections  are  to  the  line  intensity 
and  to  the  halfwidth  at  half-max  via 


S(T)  =  S(To) 


Q(To) 

Q(T) 


C2E‘(T-To) 


e 


TTo 


I  l.e-C2V,|/To  j 


and 

)(T)  =  7(To)(^)" 


(12) 

(13) 


The  values  of  Q(T)  are  obtained  from  the  polynomial  coefficients  given  in 
BLOCK  DATA  and  the  values  of  n  are  given  on  HITRAN.  For  transitions 
where  the  TIPS  are  not  yet  available  or  the  lower  state  energies  are  not  on 
the  database  (unassigncd  lines)  the  lines  are  not  written  to  the  output  file 
and  a  warning  is  issued  to  the  user  with  some  statistics  about  these  lines. 

The  program  SELFCT9I  has  been  thoroughly  tested  on  the  new 
database.  The  program  can  function  with  the  1986  database  as  well  by 
defaulting  to  the  general  search  and  using  the  subroutine  FASREAD  to 
advance  the  database.  This  does  make  selections  much  slower.  We  ha\c 
tested  all  the  options  and  enhanced  the  user  friendly  aspects  especial. y  for 
avoiding  errors.  This  version  of  vSELECT91  is  distributed  as  the  first  file  on 


3.2  Program  TR-MOM 


The  1986  HITRAN  database  saw  the  addition  of  the  transition 
moment  squared  to  the  list  of  parameters  for  each  spectral  transition.  At 
the  time  of  the  release  of  the  1986  edition  we  did  not  have  a  sufficient 
number  of  partition  functions  for  all  the  species  so  the  value  used  for  the 
transition  moment  squared  did  not  have  the  partition  sums  taken  out. 
Since  that  time  we  have  generated  the  partition  sums  for  most  of  the 
species  on  the  database  and  we  are  now  able  to  compute  the  transition 
moment  squared  from  the  line  intensity  for  most  of  the  transitions  on  the 
database.  In  addition  we  have  slightly  modified  the  parameter  on  the 
database  from  the  standard  definition  of  the  transition  moment  squared. 
This  is  to  facilitate  the  use  of  the  parameter  for  computation  of  Einstein 
coefficients.  Thus,  the  1991  edition  of  the  HITRAN  database  has  the  field 
for  the  weighted  transition  moment  squared  implemented  whenever 
possible. 

The  program  that  does  this  is  TR-MOM.  This  reads  a  line  in  the  1986 
HITRAN  format,  extracts  the  necessary  information  from  the  line,  looks  up 
the  partition  sums  and  degeneracy  factors  and  calculates  the  weighted 
transition  moment  squared,  RT^Tj',given  by 


d”  ^ 


871^ 


S(Ti)(n’)(T)  Q(T) 


(14) 


The  line  is  then  written  to  an  output  file  with  the  The  relationships 

between  Rt]ti’  and  the  Einstein  coefficients  is  described  in  section  6. 


3.3  Program  BANDSUM91 


We  have  written  a  new  version  of  the  bandsum  program  that  fully 
makes  use  of  the  vibrational  vector  coding.  This  has  allowed  us  to 
construct  a  very  general  program  that  can  be  run  cn  the  entire  database  of 
on  individual  files  by  simply  changing  the  band  center  input  selection  list. 
The  advantage  here  is  that  a  version  with  smaller  array  dimensions  can  be 
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run  on  the  PC  for  several  molecules  at  a  time.  To  run  the  program  one 
must  have  a  file  of  bands  and  band  centers.  From  this  list,  vibrational 
vectors  are  constructed  and  stored.  This  list  is  compared  to  the  vibrational 
state  of  each  line  on  the  database  as  a  function  of  molecule  and  isotope  so 
that  the  matches  can  be  stored  for  each  vibrational  band  of  each  isotopic 
species.  The  new  bandsum  program  also  does  range  checks  for  all 
parameters  of  the  database  thus  allowing  a  greater  checking  of  the 
consistency  of  the  database. 

3.4  Program  TIPS 

Program  TIPS  is  a  user  friendly  program  that  request  input  for 
Molecule,  Isotope,  and  Temperature.  The  program  uses  the  polynomial 
coefficients  for  the  fits  to  the  partition  functions,  equation  (4)  in  two 
temperature  ranges  70-400K  and  400-2005K,  or  the  coefficients  of 
equations  (5)  are  used  when  the  temperature  is  between  2005  and  300K. 
The  program  is  designed  to  use  BD-QT.for,  BD-MOL.for,  BD-ISO.for,  and  can 
easily  be  adapted  into  the  codes  of  other  users  to  determine  the  partition 
functions. 

3.5  AGENDA 

The  data  base  program  to  maintain  the  agenda  was  improved  to 
produce  the  next  version  of  HITRAN.  AGENDA  is  a  dBFile  of  all  changes  to 
be  made  to  HITRAN.  The  file  keeps  track  of  the  status  and  location  of  the 
new  data. 

The  programs  were  modified  to  make  the  screens  more  clear  for  the 
previous  options,  ADD  A  RECORD  or  MODIFY  EXISTING  RECORD.  In  addition 
two  new  options  were  added  to  CUMENU.RRG,  these  are  an  option  to 
DELETE  RECORDS  and  VIEW  RECORDS. 

In  the  ADD  A  RECORD  option  we  made  the  screen  prompt  the  user  for 
a  molecule.  If  the  entered  molecule  is  not  a  species  present  on  HITRAN  a 
message  is  issued  and  the  Enter  molecule  prompt  issued  again.  When  a 
correct  molecule  is  entered  a  routine  relates  molecule  name  to  the  number 
in  the  HITRAN  list.  In  this  fashion  the  molecule  name  and  number  always 
agree  in  the  new  database.  In  all  subsequent  options  the  molecule  name 
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and  number  can  not  be  edited  and  thus  remain  consistent.  The  fields 
Frequency,  Strength,  HW-air,  HW-self,  Energy,  Tdep-HW,  Shift,  and  Use 
have  all  been  changed  to  logicals.  The  status  line  has  been  increased  in 
size.  Also  the  record  number  of  the  record  is  displayed. 

When  modifying  records  after  each  modification  control  is  no  longer 
returned  to  the  main  menu,  one  enters  the  number  of  the  next  record  to  be 
modi^ed.  A  negative  value  will  return  control  to  the  main  menu. 

The  DELETE  RECORD  option  request  a  record  number  and  the  record 
is  then  displayed.  The  user  is  asked  if  this  is  the  correct  record,  a  Y  will 
delete  the  record  a  N  returns  the  user  to  the  main  menu.  Records  that 
have  been  deleted  can  still  be  viewed  and  edited.  In  the  upper  right  hand 
corner  of  these  records  a  message  is  given  "Record  is  marked  for 
DELETION".  VIEW  RECORDS  simply  displays  record  after  record.  A 
question.  View  next  record?  ENTER  (Y/N)  is  displayed.  The  value  Y  steps 
ahead  and  the  value  N  returns  one  to  the  main  menu. 

The  new  version  of  the  AGENDA  routine  allows  us  to  track  changes  to 
be  made  to  HITRAN  much  more  reliably. 

3.6  Program  MK-VIB 

In  preparing  to  produce  the  next  version  of  HITRAN  we  tried  to 

automate  the  procedure  as  much  as  possible.  There  are  two  main  tasks, 
first  to  remove  data  that  is  to  be  replaced  from  the  database  and  next  to 
add  the  new  and  replacement  data  to  the  remaining  lines.  The  difficulty  is 
in  designating  the  data  to  be  removed.  The  unique  specifier  is  the 
quantum  numbers  of  the  vibrational  transition  (entire  bands  are  replaced 
not  lines  within  a  band)  and  the  molecule  and  isotope  designation.  We 
have  written  a  program  that  reads  the  new  data  and  extracts  all  unique 
vibrational  labels  and  automatically  constructs  a  vibrational  table  data 

statement  for  each  species.  This  table  will  then  be  used  in  a  program  to 

remove  data  from  the  database.  The  resulting  line  list  will  contain  all 

transitions  that  are  to  be  retained  in  the  1990  version  of  HITRAN.  In 
order  to  do  this  all  data  of  each  molecule  must  be  grouped  together  and 
the  program  MK-VIB  run  on  each  molecule.  This  generates  the  list  of  all 
vibrational  transitions  of  each  isotopic  species  of  each  molecule  that  are  to 
be  removed  from  HITRAN  86. 
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3.7 


Program  REMOVE90 


This  list  of  data  to  be  removed  is  used  in  a  program,  REMOVE90,  that 
will  be  run  on  HITRAN  86.  For  each  species  of  each  molecule  data  that  will 
be  updated  is  removed  from  the  database.  This  program  counts  the 
number  of  lines  removed  and  the  number  of  lines  retained  for  each 
isotopic  species  of  each  molecule.  This  allows  a  firm  understanding  of  the 
changes  made.  The  resulting  line  file  contains  all  the  lines  on  HITRAN  86 
that  are  to  be  retained  in  the  1990  version.  This  file  will  then  be  merged 
with  the  data  files  of  new  and  updated  data  to  produce  HITRAN  1990. 

3.8  Program  ABUN 

A  program  was  written  to  take  the  most  recent^O  atomic  isotopic 
abundances  and  calculate  the  molecular  abundances  for  HITRAN.  The 
program  creates  a  table  of  abundances  and  also  prepares  a  data  statement 
for  BD-ABUN.for.  As  new  atomic  abundances  become  available  they  can  be 
added  to  ABUN  and  the  program  rerun  to  automatically  update  data  files. 

3.9  Program  Tcorr 

This  program  allows  one  to  change  the  line  intensity  to  correspond  to 
another  temperature.  It  reads  a  HITRAN  format  data  file,  and  the  user 
inputs  the  desired  temperature  of  the  output  file. 

3.10  Program  MK-REF 

This  is  a  program  that  uses  the  isotopic  species  vector  coding  to  store 
and  list  the  reference  codes  from  the  database.  This  allows  us  {o  keep 
track  of  what  reference  numbers  are  used  for  each  species  on  the 
database. 
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3.11  Module  BD-MOL 


The  molecule  chemical  symbols  are  stored  in  a  BLOCK  DATA 
subroutine  in  array  positions  0  to  32.  Any  and  all  changes  to  the  list  need 
only  be  done  here. 

3.12  Module  BD-VIB 

This  contains  several  subroutines  to  initialize  the  vibrational  vectors, 
and  to  go  back  and  forth  between  spectroscopic  notation  and  vector 
notation.  The  BLOCK  DATA  statement  contains  all  spectroscopic  vibrational 
codes  used  for  all  of  the  HITRAN  species. 

3.13  Module  BD-ISO 

This  contains  several  subroutines  to  initialize  the  isotopic  species 
vectors,  and  to  go  back  and  forth  between  the  AFGL  isotope  code  and  the 
vector  notation.  The  BLOCK  DATA  statement  contains  all  AFGL  isotopic 
species  codes  used  for  all  of  the  HITRAN  species. 

3.14  Module  BD-ABUN 

This  contains  a  data  statement  in  a  BLOCK  DATA  routine  for  the 
isotopic  abundance  of  each  molecular  species  on  the  database. 

3.15  Module  BD-QT 

This  is  a  BLOCK  DATA  subroutine  that  contains  the  polynomial 
coefficients  for  the  three  temperature  ranges  the  partition  sums  were 
calculated  for  70-400K,  400-2005K,  and  2005-3000K.  The  partition  sum 
for  each  isotopic  species  is  stored  at  the  reference  temperature  296K.  And 
the  state  dependent  and  independent  degeneracy  factors  are  stored. 
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4.0  1991  HITRAN 


The  work  on  the  new  data  to  the  1991  HITRAN  is  briefly 
summarized  below.  Some  new  species  have  been  added  to  the  HITRAN 
database.  To  keep  track  of  these  AGENDA  has  been  modified  to  consider 
these.  The  new  species  are  COF2  (molecule  #  29),  SF6  (molecule  #  30),  H2S 
(molecule  #  31),  and  HCOOH  (molecule  #  32).  This  allows  us  to  track 
additions  and  updates  for  these  species. 

A  fifth  isotopic  species  for  carbon  monoxide  was  added  to  the  1991 
database,  This  has  shifted  the  isotopic  vector  positioning.  All 

programs  using  the  isotopic  vector  structure  and  the  BD-ISO  routines  were 
modified  to  account  for  this.  The  isotopic  abundance  program  was 
modified  and  the  new  abundances  created.  These  are  given  in  Table  2  with 
the  appropriate  number  of  significant  figures. 

4. 1  H2O 

Water  vapor  has  seen  some  significant  improvements  for  this  edition 
of  the  database.  We  have  added  the  data  of  Flaud  et  al.5i-53  for  the 
regions  8036-9482  cm'^,  9603-1 1481cm-l,  and  13238-22657  cm'^,  the 
data  of  Toth54.55  vvas  used  for  the  5904-12741  cm'l  region  and  for  the  v2 
region.  The  line  count  for  these  files  is  1874,  2484,  4608,  3637,  and  4618 
lines  respectively.  For  both  of  the  Toth  files  strong  lines  were  not 
included  so  they  had  to  be  taken  from  the  1986  database.  This  is 
discussed  below.  In  general  for  this  data  we  had  to  perform  the  following 
task 

1 .  Reformat  the  data  to  the  H1TRAN86  format. 

2.  Convert  line  intensities  from  cm'^  atm*^  units  to  cm'Vmolec  cm'2. 

3.  Convert  line  intensities  to  300°  K  to  reference  temp  (296°  K). 

4.  Add  the  intensity  error  codes. 

5.  Add  the  lower  state  energies  for  each  line. 

6.  Add  reference  codes. 

7.  Add  halfwidths  and  temperature  exponents 

8.  Add  lower  state  energies. 
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Table  2.  Isotopic  Abundances  for  Species  on  HITRAN91 


Molecule 

HITRAN 

isotope 

code 

Abundance 

Molecule 

HITRAN 

isotope 

code 

Abundanc 

H20 

161 

.997317 

HF 

19 

.99984425 

181 

.00199983 

171 

.000372 

HCI 

15 

.757587 

162 

.00031069 

17 

.242257 

002 

626 

.98420 

HBr 

19 

.50678 

636 

.01106 

11 

.49306 

628 

.0039471 

627 

.000734 

HI 

17 

.99984425 

638 

.00004434 

637 

.00000825 

00 

56 

.75591 

828 

.0000039573 

76 

.24172 

728 

.00000147 

OCS 

622 

.93739 

03 

666 

.992901 

624 

.04158 

668 

.00398194 

632 

.01053 

686 

.00199097 

822 

.0018797 

N20 

446 

.990333 

H2CD 

126 

.98624 

456 

.0036409 

136 

.01108 

546 

.0036409 

128 

.0019776 

448 

.00198582 

447 

.000369 

HOO 

165 

.75579 

167 

.24168 

CD 

26 

.98654 

36 

.01108 

N2 

44 

.9926874 

28 

.0019782 

27 

.000368 

HCN 

124 

.98511 

38 

.00002222 

134 

.01107 

125 

.0036217 

CH4 

211 

.98827 

311 

.OHIO 

CH3CI 

215 

.74894 

212 

.00061575 

217 

.23949 

O2 

66 

.995262 

H2O2 

1661 

.994952 

68 

.00399141 

67 

.000742 

C2H2 

1221 

,97760 

1231 

.02197 

NO 

46 

.993974 

56 

.0036543 

C2H6 

1221 

.97699 

48 

.00199312 

PH3 

nil 

.99953283 

SO2 

626 

.94568 

646 

.04195 

COF2 

269 

.98654 
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Table  2.  continued 


N02 

646 

.991616 

SF6 

29 

.95018 

NH3 

4111 

.9958715 

H2S 

121 

.94988 

Sill 

.0036613 

HOCX)H 

126 

.98390 

HN03 

146 

.989110 

OH 

61 

.997473 

81 

.00200014 

_ 62 _ 

Dr.  Brown56  informed  us  that  the  near-IR  and  visible  data  have 
many  of  the  strong  lines  missing  due  to  saturation.  The  same  was  found  in 
the  V2  region  data  from  Toth.  To  alleviate  this  problem  data  for  the 
following  bands  were  extracted  from  the  1986  HITRAN,  (OIO)-(OOO), 
(040)-(000),  (120)-(000),  r021)-(000),  (200)-(000),  (lOl)-(OOO),  (002)- 
(000),  (131)-(000),  (310)-(000),  (211)-(000),  (112)-(000),  and  (013)- 
(000).  This  data  were  compared  to  that  of  Toth  and  all  lines  not  present  in 
the  Toth  data  were  retained.  This  gave  us  2307  lines  in  the  near-IR 
region,  669  lines  in  the  visible  region,  and  86  lines  for  the  V2  region  that 
will  be  added  to  the  1991  database  to  complete  the  data  for  these  bands. 

4.2  CX>i 

Dr.  Laurence  Rothman  created  a  file  containing  all  CO2  lines.  This  file 
has  the  most  up  to  date  halfwidths,  and  temperature  exponents,  error 
codes,  references  and  transition  moments. 

4.3  O3 

Permission  was  granted  by  Academic  Press  to  use  the  ozone  data  of 
Flaud  et  al.57  on  the  1990  version  of  HITRAN.  Curt  Rinsland  of  NASA 
Langley  sent  us  the  complete  data  for  ozone  in  six  compressed  files.  The 
format  of  the  files  is  a  modified  1982  format.  I  wrote  a  program  to 
reformat  the  data  into  the  1986  format,  add  the  self  and  air-broadened 
halfwidths,  the  temperature  dependence  of  y,  add  the  error  code  for  y,  and 
to  add  data  to  the  reference  code  as  requested  by  Academic  Press. 
Rinsland  suggested  not  adding  position  and  intensity  error  codes  at  this 
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time.  The  pure  rotation  data  has  only  4  significant  figures  after  the 
decimal  whereas  we  allow  for  6.  He  suggest  using  the  new  file  due  to  the 
improvements  in  the  intensities  and  commented  that  the  four  significant 
figures  was  already  more  than  the  data  could  support. 

The  reformatting  program  was  run  on  all  the  files.  The  reformatted 
files  were  merged  in  one  file  called  OZONE.90  which  contained  168881 
lines.  Inspection  of  the  data  statement  revealed  the  ozone  data  would  be  a 
complete  replacement  of  ozone  data  on  HITRAN  86.  The  new  data  now 
contains  66  bands  for  the  principal  species  and  3  bands  each  for  the  lesser 
j;.otopic  species. 

4.4  CD 

The  carbon  monoxide  files  were  obtained  from  R.  Tipping^S.  There 
are  5  files,  one  for  each  isotopic  species.  Error  codes  and  reference  codes 
were  added  to  the  files.  The  data  are  in  the  files  C051.91,  C052,91, 
C053.91,  C054.91,  and  C055.91  and  they  contain  800,  700,  700,  700,  and 
70*;  lines  respectively. 

4.5  CH4 

All  of  the  methane  on  the  database  was  replaced  with  a  file  sent  to 
us  by  Linda  Brown^^. 

4.6  NO 

John  Ballard’s  measurements^^  of  the  NO  fundamental  will  be  added 
to  the  database.  Ballard  has  12"^  from  P  13.5  to  R  19.5.  The 

fundamental  was  extracted  from  the  1986  database  to  get  the  line  missing 
in  Ballard’s  data.  Of  the  417  lines  from  HITRAN  1986,  there  were  289 
lines  retained.  The  error  and  reference  codes  were  set  to  zeros  for  these 
lines. 
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4.7  HNOs 


There  are  six  files  for  HNO3  from  Aaron  Goldman^O.  These  lines  are 
in  the  HITRAN  format.  The  only  change  needed  was  the  addition  of  the 
references.  The  data  are  in  the  files  HN03N2.91,  HN03N6.91,  HN03N7.91, 
HN03N8.91,  HN03N9.91,  and  HN03N89.91  and  they  contain  32340,  8379, 
8013,  7101,  8091,  and  9709  lines  respectively. 

4.8  Hydrogen  Halides 

Six  files  were  obtained  from  R.  Tipping^^  for  the  hydrogen  halides, 
one  each  for  HF  and  HI  and  two  each  for  HCl  and  HBr.  The  references  and 
error  codes  were  added.  The  HCl  files  Tipping  sent  were  in  WORDSTAR 
format  and  had  to  be  converted  to  text.  One  of  these  files  had  a  decimal 
point  in  the  wrong  field  indicating  that  all  or  part  of  the  data  was  entered 
by  hand(This  may  haunt  us  in  the  future).  The  final  files  for  this  data  are 
HF.91,  HC1151.91,  HC1152.91,  HBr61.91,  HBr62.91,  and  HI.91  and  they 
contain  107,  203,  168,  200,  198,  and  237  lines  respectively. 

4.9  N2 

For  nitrogei:  the  data  from  HITRAN  86  was  extracted  and  the 
intensities  scaled  by  1.049  and  the  halfwidth  set  to  0.047  cm' l /atm  as 
recommended  by  Curt  Rinsland^*.  Reference  codes  were  added  to  make 
the  final  file  of  117  lines  called  N2.91. 

4.10  H2O2 

The  V6  band  was  obtained  from  John  Hillman^2  to  replace  the  data  on 
HITRAN  86.  The  file  was  reformatted  and  the  shift  and  self-broadened 
halfwidth  were  set  to  0,  Vair  set  to  0.10  cm'* /atm,  n=0.5,  and  reference 
codes  added.  The  final  file  is  H2O2.9I  and  contains  4561  lines. 
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4.11  C2H2 


We  have  119  lines  of  data  for  the  V5  band  of  C2H2^^  this  is  new  to 
the  database.  The  error  codes  appeared  to  be  in  the  style  of  the  old  date 
code  (288)  meaning  February  1988.  This  was  changed  to  000.  Reference 
codes  were  added  and  the  final  file  is  C2H2.9I. 

4.12  C2H6 

This  is  the  V7  band  of  C2H6^^,  it  is  new  data  to  HITRAN.  There  are 
421  lines  of  data.  Reference  codes  were  added  to  the  data.  The  data  are  in 
the  file  C2H6.9I. 

4.13  COF2 

This  is  a  new  molecule  to  the  database.  It  has  been  assigned  as 
molecule  29  and  the  isotope  code  is  1  corresponding  to  the  GL  code  of  269, 
shorthand  for  l2cl6o^^F2.  There  are  three  files  for  the  V2,  V4  ,  and  V6  bands 
from  Goldman  et  al.^^  and  they  contain  7234,  3250,  and  7758  lines 
respectively. 

4.14  SF6 

This  is  a  new  molecule  to  the  database.  It  has  been  assigned  as 
module  30  and  the  isotope  dode  1  corresponds  to  the  GL  code  29,  short 
hand  for  32s '^F6.  The  data  is  from  L.  Brown^^  and  is  for  the  V3  band. 
These  are  a  total  of  11520  lines  for  this  data. 

4.15  Reference  Codes 

Reference  codes  have  been  implemented  for  all  the  new  and  updated 
data  for  the  1991  edition  of  HITRAN.  To  label  data  from  the  previous 

version  the  error  and  reference  codes  have  all  been  set  to  zeros.  This  will 
indicate  the  reference  of  Rothman  et  al.L  A  file  of  the  references  will  be 
one  of  the  files  distributed  with  the  HITRAN  database.  The  1991  edition  of 
HITRAN  will  have  both  the  error  and  reference  codes  fully  active  for  the 
first  time. 


29 


4.16  1986  HITRAN  Data  for  the  1991  Edition 

From  a  careful  inspection  of  the  AGENDA  output  listing  we  have 
determined  which  molecular  species  will  be  updated  in  the  1991  edition  of 
HITRAN.  A  procedure  has  been  developed  to  remove  from  the  1986 
database  all  the  data  corresponding  to  the  new  data.  This  will  leave 
behind  all  the  1986  data  that  will  be  retained  in  the  1991  edition.  For  the 
molecules  CO2,  O3,  CO,  CH4,  and  the  hydrogen  halides  all  data  is  to  be 
replaced  thus  these  lines  are  simply  removed  from  the  1986  database.  For 
the  other  species  to  be  updated,  H2O,  NO,  HNO3,  N2,  H2O2,  C2H2,  and  C2H6  the 
following  procedure  was  used.  For  each  species  all  data  was  grouped  into 
one  file.  The  file  was  then  run  through  MK-VIB  which  finds  all  unique 
vibrational  bands  in  the  data  and  creates  a  data  statement  of  the 
vectorized  transitions  for  each  isotopic  species,  e.g.  for  water  vapor  the  V2 
transition  is  given  by  (OIO)-(OOO),  the  vector  code  is  2  to  1  and  the 
vectorized  transition  is  given  by  1000*IVUP+IVLO  thus  for  V2  we  have 

2001.  This  generates  a  data  statement  of  all  vibrational  bands  for  each 
isotopic  species  of  a  molecule  that  must  be  removed  from  the  1986 

database.  Nothing  must  be  typed  by  hand  in  this  procedure.  The  data 

statements  for  each  species  were  then  put  into  a  program  REMOVE91 

which  was  run  on  the  1986  database.  The  line  counts  from  this  procedure 
are  given  in  Table  3  and  they  agree  exactly  with  what  one  would  expect 
from  AGENDA. 


Table  3  Line  Counts  from  REMOVE91 


MoJ _ ISQ...  rciaingd _ r.gmo..ygd _ Mol  iso  retained  removed 


H20 

161 

12098 

17268 

HNO3 

146 

28496 

7492 

181 

5076 

1039 

171 

2662 

754 

OH 

61 

8521 

0 

162 

8305 

0 

81 

65 

0 

62 

90 

0 

C02 

626 

0 

26201 

636 

0 

8766 

HF 

19 

0 

62 

628 

0 

13046 

627 

0 

6624 

HCl 

15 

0 

103 

638 

0 

2322 

17 

0 

97 

637 

0 

1585 

828 

0 

722 

HBr 

19 

0 

128 

728 

0 

288 

11 

0 

128 

03 

666 

0 

47938 

HI 

17 

0 

145 

668 

0 

1149 

686 

0 

993 

00 

56 

3011 

0 

76 

.3009 

0 

N20 

446 

17792 

0 

456 

1876 

0 

OCS 

622 

461 

0 

546 

1901 

0 

624 

99 

0 

448 

1774 

0 

632 

93 

0 

447 

782 

0 

822 

84 

0 

(D 

26 

0 

266 

H2OO 

126 

1772 

0 

36 

0 

137 

136 

563 

0 

28 

0 

111 

128 

367 

0 

27 

0 

60 

HOO 

165 

8057 

0 

CH4 

211 

0 

14587 

167 

7508 

0 

311 

0 

1130 

212 

0 

2057 

N2 

44 

0 

117 

Tabic  3 

continued 

02 

66 

1132 

0 

HCN 

124 

703 

0 

68 

612 

0 

134 

34 

0 

67 

510 

0 

125 

35 

0 

NO 

5590 

417 

CHja 

215 

3977 

0 

56 

699 

0 

217 

2710 

0 

48 

679 

H2O2 

1661 

883 

2389 

SO2 

626 

23372 

0 

646 

287 

0 

C2H2 

1221 

1052 

0 

1231 

87 

0 

NO2 

646 

20067 

0 

C2H6 

1221 

4328 

0 

NH3 

4111 

4856 

0 

■S.in _ 

%\  - 

_ 0 _ 

pmnii _ 

-2886 _ 

_ Q 

4.17  Structure  of  Cross-Section  Files  on  HITRAN 

The  structure  of  the  cross-sectional  files  on  HITRAN  had  never  been 
fixed  to  a  precise  format.  For  the  1991  version  >ve  have  decided  to  use  the 
following  format;  each  file  will  start  with  a  100  byte  header  which  contains 
the  chemical  symbol  right  justified  in  AlO  format,  the  initial 
wavenumber  and  final  wavenumber  each  in  flO.4  format,  the 
number  of  points  of  data  in  110  ,  the  temperature  (in  Kelvins)  and 
pressure  (in  atmospheres)  corresponding  to  the  data  in  FI 0.4,  the 
Maximum  cross-section  in  the  data  in  E10.3,  and  source  of  the  data 
with  comments  in  A30.  The  cross-sections  corresponding  to  this  header 
follows  immediately  in  100  byte  lines  (10E10.3)  at  a  wavenumber  spacing 
defined  by 

Av  =  (v final  ■  Vinitial)/Npoinis  (15) 
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The  cross-sections  and  the  Xmax  are  in  units  of  cm'V(niolec*cin*2)  at  the 
temperature  T.  If  the  data  ends  in  the  middle  of  a  line,  the  line  is  zero 
filled  to  the  end.  See  Figure  6  for  details.  The  next  group  of  cross-.section 
data  begins  on  the  next  line  with  the  header.  This  file  will  be  put  on 

HITRAN  91  in  blocked  form,  5100  bytes/block.  This  means  that  many  of 

the  cross-section  records  will  begin  within  the  block  structure  as  shown  in 
Figure  6. 

4.18  The  1991  edition  of  the  HITRAN  Database 

It  should  be  noted  that  several  iterations  took  place  to  produce  the 

1991  edition  of  the  HITRAN  database.  Errors  that  were  detected  in  the 

data  were  corrected.  This  included  among  other  corrections;  the  HCl  data 
from  Tipping,  addition  of  other  data  for  NO2,  OCS,  SF6,  and  HNO3,  resetting 
reference  codes  for  the  Camy-Peyret  et  al.  visible  region  H2O  data,  NO2 
duplicate  lines,  etc.  Each  time  an  iteration  was  done  we  had  to  recreate  the 
blocked  database  with  the  lookup  vectors.  This  was  tedious,  required 
running  several  programs,  and  files  had  to  be  created  and  deleted  during 
the  process.  To  make  the  task  more  efficient,  the  programs  were 
consolidated  into  one  called  MAKE_BLOCKED_HITRAN.for.  To  run  this,  one 
first  must  sort  the  unblocked  HITRAN  line  file.  This  must  be  done  by 
defining  scratch  working  files  by  issuing  the  following  commands  on  the 
VAX, 

$ASSIGN  SCRATCH  SORTWORKO 
SASSIGN  SCRATCH  SORTWORKl 


SASSIGN  SCRATCH  SORTWORK9 

SORT/STABLE/KEY=(POS:4,S1ZE:12)AVORK_F1LES=10  filein  filcoiit. 
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Header  •  •  •  N  points  •  •  • 


data 

51x100  bytes 


data 

51x100  bytes 


Block  1 


Block  2 


Block  J 


^  Next  group  of 
X-section  data 


data 

51x100  bytes 


Block  J+1 


Figure  6.  Structure  of  the  1991  Cross-section  Files. 


MAKE_BLOCKED_HITRAN  then  takes  the  sorted  output  file  and  blocks  it 
into  a  scratch  filefHlTRAN  format),  the  file  is  rewound  and  the  frequency- 
block  table  is  produced.  From  this  file  the  frequency-block 
correspondence  table  is  created,  currently  this  puts  26  blocks  per  marker, 
and  the  edition  number  is  entered  and  the  correspondence  blocks  created 
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(there  are  two  such  blocks).  Finally,  the  correspondence  blocks  and  the 
blocked  HITRAN  file  are  written  to  one  file,  the  new  edition  of  the  HITRAN 
database.  The  final  database  we  created  was  labeled  HITRAN  Edition  03 
1991  and  it  contains  647672  lines  in  total.  In  addition  there  are  102  = 
2x51  "lines"  of  header  and  position  pairs.  Table  4  gives  the  line  count  for 
each  isotopic  species  on  the  database. 

New  on  this  database  is  the  use  of  the  reference  fields  for  the 
position,  intensity,  and  air-broadening  coefficient.  We  have  made  full  use 
of  these  fields  for  all  the  new  data  and  have  set  the  remaining  1986 
database  lines  to  have  a  reference  code  of  0  for  each  parameter.  In  the 
new  data  we  have  also  tried  to  implement  the  error  codes  whenever 
possible.  The  error  codes  for  the  1991  database  are  given  in  Table  5  for 
the  oosition  and  in  Table  6  for  the  line  intensity  and  the  haltwidth. 

5 . 0  BROADENING  COEFHCIENTS 

Much  work  was  done  in  the  area  of  the  broadening  coefficients  for 
the  1991  database.  This  involved  both  calculations  and  better  use  of  the 
experimentally  measured  values  for  several  of  the  species  on  the  database. 

5.1  Air-Broadened  Halfwidths  for  Water  Vapor  on  HITRAN 

Dr.  Curtis  Rinsland  of  NASA  Langley  sent  a  preprint's  on 
measurements  of  air-broadened  halfwidths  of  water  vapor  he  made  with 
Aaron  Goldman,  Mary  Ann  Smith,  and  Malathy  Devi.  I  reviewed  the 
preprint  and  was  surprised  by  the  agreement  (or  lack  of  agreement) 
between  the  HITRAN  values  and  their  measurements.  I  knew  there  would 
be  differences  but  I  did  not  expect  them  to  be  so  large.  When  the 
calculations  were  done,  1979-80,  there  were  not  many  measurements  to 
compare  with.  For  example,  comparison  made  at  that  time  for  J<8  pure 
rotation  lines  gave  an  average  absolute  percent  difference  of  12.2  for  some 
14  lines.  I  think  the  differences  found  in  the  paper  of  Rinsland  et  al.  are 
understandable  given  the  nature  of  the  work  done  in  1980. 
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Table  4  Line  Count  for  HITRAN91 


[HHli 

H20 

161 

48523 

181 

6357 

171 

3744 

162 

8305 

002 

626 

26887 

636 

8892 

628 

13343 

627 

6625 

638 

2312 

637 

1584 

828 

721 

728 

288 

03 

666 

142221 

668 

19147 

686 

7513 

N20 

446 

17792 

456 

1876 

546 

1901 

448 

1774 

447 

782 

(D 

26 

800 

36 

700 

28 

700 

27 

700 

38 

700 

CH4 

211 

35588 

311 

4926 

212 

6457 

O2 

66 

1132 

68 

612 

67 

510 

NO 

46 

6007 

56 

699 

48 

679 

SO2 

626 

23372 

646 

287 

NO2 

646 

26296 

NH3 

4111 

4856 

5111 

961 

Iso 

#  transitions 

HNO3 

146 

143021 

OH 

61 

8521 

81 

65 

62 

90 

HF 

19 

107 

HCl 

15 

203 

17 

168 

HBr 

19 

200 

11 

198 

HI 

17 

237 

(30 

56 

3011 

76 

3009 

0(3 

622 

461 

624 

99 

632 

93 

822 

84 

H2OO 

126 

1772 

136 

563 

128 

367 

Hoa 

165 

8057 

167 

7508 

N2 

44 

117 

HCN 

124 

703 

134 

34 

125 

35 

CH3CI 

215 

3977 

217 

2710 

H2O2 

1661 

5444 

C2H2 

1221 

1 172 

1231 

86 

C2H6 

1221 

4749 

PH3 

nil 

2886 

COF2 

269 

18242 

. 29 

_ 11520 
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Table  5.  Error  Code  for  the  Wavenumber 


Code 

Error  Range 

0 

>  1./  or  undefined 

1 

>  .1  and  <  1. 

2 

>  .01  and  <  .1 

3 

^  .001  and  <  .01 

4 

>  .0001  and  <  .001 

5 

>  .00001  and  <  .0001 

6 

<  .00001 

Table  6.  Error  Code  for  the  Line  Intensity  and  Halfwidth 


Code 

Error  Range 

0 

Undefined 

1 

Default  or  Constant 

2 

Average  or  Estimate 

3 

^  .20% 

4 

^  10%  and  <  20% 

5 

^  5%  and  <  10% 

6 

>  2%  and  <  5% 

7 

>  1%  and  <  2% 

8 

<1% 

The  values  presented  on  HITRAN  are  from  two  sources.  One  is  the 
database  of  N2 -broadened  halfwidths  presented  in  the  Applied  Optics 
article^^  and  the  others  are  from  an  averaging  algorithm  devised  by  Dick 
Davies^^.  The  calculations  followed  the  work  of  reference  68.  These 
calculations  only  considered  the  dipole  (H2O)  -  quadrupole  (N2)  interaction 
and  hence  will  be  too  small  because  the  quadrupole  (H2O)  -  quadrupole 
(N2)  interaction  has  been  omitted.  In  Davies'  work  the  QFT  scaling 
parameter  a  was  set  to  2.79  to  match  a  calibration  line  (5  2  3->6  1  6)  and 
the  parameter  bmin  set  to  1.5  A.  We  have  since  found  a  =  2.50  to  be  a 
better  choice.  At  the  time,  Davies  was  trying  to  adjust  the  theory  to 
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explain  some  very  narrow  line  widths  for  high  J  lines  that  were  measured. 
This  was  accomplished  by  the  low  value  of  bmin-  As  stated  by  Davies  and 
Oli^S,  this  choice  of  bmin  give  questionable  results  for  intermediate  J 
values.  This  is  what  is  observed.  Note  Benedict  and  Calfee^^  used  a  bmin 
value  of  3.2  A  thus  there  results  will  be  proportionally  larger.  It  appears 

that  the  optimization  of  the  calculations  for  the  high  J  lines  has  led  to 

halfwidths  that  are  consistently  low  as  shown  by  the  results  of  Rinsland  et 
al..  The  other  quantities  that  enter  into  the  calculation,  energies,  reduced 
dipole  and  quadrupole  line  strengths,  etc.  were  also  not  as  well  known  for 

the  intermediate  to  upper  states  of  water  vapor  leading  to  more 

uncertainty  in  these  calculations. 

We  did  a  few  calculations  to  compare  with  Tables  I  and  IV  of 
reference  65.  The  calculations  are  essentially  the  method  used  in  the  1988 
paper  on  temperature  dependence  by  Gamache  and  Rothman^O  and  is 
nearly  equivalent  to  that  of  Tejwani^l.  In  ref.  70  the  Tipping  tr.\,ectory 
model^^  vvas  used  and  Tejwani  uses  that  of  Robert  and  Bonamy^3  |t  was 
shown  for  ozone^'*  that  these  difrer  by  3%  with  the  Tipping  model  giving 
the  larger  values.  The  calculations  done  here  are  somewhat  larger  than 
Tejwani's  and  without  knowing  the  explicit  molecular  constants  he  used  we 
cannot  address  the  differences.  What  is  more  important  is  that  although 
there  appears  to  be  agreement  for  these  lines  for  the  other  pure  rotational 
transitions  studied  there  remains  large  differences.  For  one  of  the 
measured  lines  (6  2  5-»7  5  2)  we  had  done  a  QFT-ID  (quantum  Fourier 
Transform  Theory  with  Improved  Dynamics)  calculation  some  years  ago 
and  the  agreement  is  somewhat  better  than  the  BKG-ID  value.  Comparing 
with  the  14  experimental  lines  mentioned  above  the  QFT-ID  method  gave 
an  average  absolute  percent  difference  of  7.7.  However  as  higher  J  states 
are  considered  the  agreement  becomes  worse.  Clearly  work  in  needed  on 
comparing  and  developing  the  theories  and  the  intermolecular  potentials. 

5.1.1  Direct  Access  Files  for  H2O  and  O3  halfwidths 

For  the  updating  of  halfwidths  on  the  database  for  O3  and  H2O  a  new 
approach  was  developed.  This  is  the  use  of  direct  access  files  on  the  PCs  or 
the  VAX  to  update  the  halfwidths.  The  approach  is  similar  to  the  mass 
storage  approach  that  was  used  on  the  CYBER  however  it  is  done  on  the  PC 
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or  VAX  with  much  greater  conirol.  We  have  created  direct  access  files  of 
N2  and  self  broadened  halfwidths  of  H2O  which  will  run  in  ADHWH20. 
From  my  calculations  1  have  obtained  the  average  halfwidth  as  a  function 
of  J  which  has  also  been  added.  These  are  given  in  Table  7. 


Table  7.  Average  N2-  and  Self-Broadened  Halfwidths  of  H2O 


T" 

XfVl 

(X 

a. 

0 

.10554 

1 

.51760 

1 

1 

.11038 

7 

.53498 

7 

2 

.10532 

1  9 

.53558 

1  9 

3 

.09897 

33 

.51965 

3  3 

4 

.09217 

47 

.48387 

47 

5 

.08512 

6  1 

.44795 

61 

6 

.07762 

75 

.41578 

7  5 

7 

.07016 

89 

.36937 

89 

8 

.06319 

103 

.32831 

103 

9 

.05690 

117 

.30085 

1  17 

10 

.05139 

131 

.28184 

131 

1  1 

.04650 

145 

.26001 

145 

1  2 

.04199 

159 

.23730 

159 

1  3 

.03787 

173 

.21433 

173 

1  4 

.03405 

187 

.19057 

1  87 

15 

.03038 

201 

.17076 

201 

1  6 

.02701 

215 

.15836 

215 

1  7 

.02405 

229 

.14736 

229 

1  8 

.02277 

208 

.14874 

208 

1  9 

.02045 

226 

.13518 

226 

20 

.01856 

156 

.12327 

1  56 

For  ozone  I  took  the  database  of  theoretical  values'^-*’  and  stored  these 
in  a  direct  access  file.  The  program  to  add  the  halfwidths  operates  as 
follows:  given  an  initial  and  final  state  the  DAF  is  searched  for  data.  If 
available  the  N2  broadened  value  is  used.  When  there  is  no  calculation 


and  J"<40  the  averaged  theoretical  halfwidth  is  used.  For  these  cases,  the 
value  is  then  scaled  by  0.95  to  get  the  corresponding  air-broadened  value.. 
It  is  then  scaled  another  9%  based  on  M.  Smith's  comparisons^  of  the 

calculations  to  experiment.  For  air-broadening  above  J"=40  Curt 
Rinsland's  polynomial  values^S  are  used.  For  the  self-broadened  values 
the  polynomial  expression  of  Smith  et  al.S8  determined  by  a  fit  to  355 
measurements  is  used. 

The  addition  of  air-  and  self-broadened  halfwidths  and  temperature 

exponents  for  H2O  has  undergone  a  major  change.  The  previous  approach 
used  a  database  of  calculated  values  which  are  known  to  be  in  error 

especially  for  the  intermediate  J  lines.  The  new  approach  uses 
experimentally  determined  halfwidths  whenever  possible.  The  procedure 
is  as  follows:  for  each  perturber  a  master  file  was  created  consisting  of  the 
vibrational  and  rotational  quantum  numbers,  the  halfwidth,  the  error  in 

the  halfwidth,  the  shift,  and  a  reference  number.  These  files  are  used  as 
input  to  a  program  to  create  the  direct  access  files  (DAFs)  for  y  and  n. 
Currently  the  DAFs  of  halfwidths  contain  the  value  of  the  halfwidth,  its 
error  code  and  reference;  the  DAFs  of  temperature  exponents  contain  n 
and  the  reference..  In  creating  direct  access  files  the  experimental  data 

are  given  preference.  When  no  experimental  data  exist  scaled  calculated 

values  of  Gamache  and  Davies^^  are  used.  The  scaling  values  were 

determined  separately  for  each  broadening  gas  by  comparing  the 
theoretical  and  experimental  values.  From  this  analysis  scale  factors  as  a 
function  of  J  were  determined  for  both  perturbers.  The  master  files  of  the 
halfwidths  can  be  updated  at  any  time  and  the  new  DAF  created. 

5.2  Collision  Broadened  Halfwidths  for  CO2 

Ro.senmann  et  al.^*’-^*  have  computed  the  collision  broadened 
halfwidths  and  temperature  dependence  for  CO2  perturbed  by  CO2.  H2O.  NS. 
and  O2.  This  data  was  taken  and  used  in  a  subroutine  to  evaluate  the 

halfwidth  for  gaseous  mixtures  and  temperature  exponent  for  the 
temperature  dependence  of  the  halfwidth  for  CO2  lines.  Ihis  will  be  added 
to  the  CO2  line  generator  program  of  Rothman  to  evaluate  air-broadened 
halfwidths  and  the  temperature  exponent.  Although  it  is  a  consistent  set 
of  data  we  have  found  it  to  be  greater  than  the  experimental  results  of 
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Johns *2  and  Dana  et  al.83  for  all  values  of  m.  We  have  addressed  the 
assumption  of  replacing  the  velocity  integral  by  the  mean  relative  thermal 
velocity,  MRTV,  as  a  possible  cause  for  the  differences.  This  approximation 
becomes  more  suspect  at  elevated  temperatures.  The  reason  for  this  is 
that  at  elevated  temperatures  the  Maxwell-Boltzmann  distribution  of 
velocities  spreads  out  and  the  assumption  of  replacing  the  integral  over 
velocities  by  the  mean  relative  thermal  velocity  is  less  valid.  From  our 
calculations  we  have  found  that  the  agreement  is  better  at  elevated 
temperatures.  We  have  sought  to  explain  this  and  have  found  some 
interesting  results. 

The  Baranger-Kolb-Griem  formalism^^  gives  the  Lorentz  component 
halfwidth  at  half  maximum  (HWHM),  Vif,  of  isolated  lines  as  the  real  part  of 
the  diagonal  matrix  elements  of  the  resolvent  operator.  Within  the  impact 
approximation,  the  broadening  coefficient  is  given  by 


Yi/(cm-Vatn) 


I  2rtbSii(b.u.J2)cto 


(16) 


where  N  is  the  number  density  of  perturbers,  v  the  relative  velocity  of  the 
colliding  pair,  b  is  the  impact  parameter,  c  the  speed  of  light,  pj2  is  the 
Boltzmann  factor  for  the  perturber  state  J2»/(v)  is  the  Maxwell-Boltzmann 
velocity  distribution  function,  and  Sif(b,v,  J2)  has  been  called  the  collision 
efficiency  function  or  the  diffusion  operator.  It  is  found  that  by  replacing 
the  velocity  integral  by  the  mean  relative  thermal  velocity  and  evaluating 
the  diffusion  operator  at  P,  gives  results  generally  within  5-10%  of  that 
calculated  by  Eq(16).  The  the  broadening  coefficient  in  terms  of  the  mean 
relative  thermal  velocity,  P,  is  given  by 


Y/7{(nTyatm) 


2xbSii(b,v,J2)  cto 


(17) 


Formally  this  expression  is  developed  by  assuming  that  Sif(b.  v,  J2)  does 
not  depend  on  velocity,  thus  allowing  the  integral  on  velocity  to  be  done 
analytically.  However  the  diffusion  operator  is  known  to  vary  with 


velocity.  The  reason  equation  (17)  is  used  in  place  of  (16)  is  the  savings  in 
time  for  the  calculation. 

The  need  for  broadening  coefficients  with  an  uncertainty  of  less  than 
10%  has  prompted  us  to  study  the  approximation  leading  to  Eq.  (17).  We 
have  investigated  the  velocity  dependence  of  Sif(b,  v,  J2)  for  CO2  perturbed 
by  N2,  O2,  and  CO2  and  found  that  the  agreement  between  values  calculated 
via  equations  (16)  and  (17)  is  fortuitous  and  conclude  that  formally  the 
velocity  integral  should  be  explicitly  evaluated.  Given  this,  we  have 

recalculated  the  halfwidths  and  temperature  exponents  for  N2,  O2,  and  CO2 
broadening  of  CO2  by  explicitly  using  equation  (16). 

We  have  compared  these  values  to  the  calculations  of  Rosenmann  et. 
ai.80.81.  In  general  we  find  that  from  Iml  up  to  20  our  calculations  are 
lower  than  Rosenmann’s,  bringing  them  closer  to  the  experimental  values. 
The  reason  for  the  lowering  is  that  the  MRTV  calculations  have  an  artificial 
resonance  in  the  diffusion  operator  from  «  5  xl0*8  to  7  xlO'*  cm.  This  can 
be  seen  in  figure  7  for  the  J“=l  R-branch  transition  of  CO2  interacting  with 
the  J2=0  level  of  N2  at  296K.  At  lml=20  the  curves  cross  and  our 

calculations  are  lower  than  Rosenmann’s  by  3%  at  lml=25  decreasing  to 
0.5%  for  lml>100.  It  should  be  emphasized  that  the  MRTV  calculations 
even  when  they  agree  with  the  full  calculation  do  so  only  by  a  fortuitous 
cancellation  of  errors.  Differences  in  the  interruption  functions 

approaching  90%  (Figure  8)  have  been  observed  between  the  two  methods. 

6.0  LINE  INTENSITY.  TRANSITION  MOMENT  SQUARED.  EINSTEIN 

COEFFICIENTS.  AND  OSCILLATOR  STRENG'fi^S 

The  1986  HITRAN  database  saw  the  addition  of  the  transition 

moment  squared  to  the  list  of  parameters  for  e  ich  spectr.il  transition.  At 

the  time  of  the  release  of  the  1986  edition  we  did  not  have  a  sufficient 
number  of  partition  functions  for  all  the  species  so  the  value  used  for  the 
transition  moment  squared  did  not  have  the  partition  sums  taken  out. 
Since  that  time  we  have  generated  the  partiti  m  sums  for  most  of  the 
species  on  the  database  and  we  are  now  able  to  compute  the  transition 
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Figure  7.  Comparison  of  diffusion  operator  calculated  from  equations 
(16)  and  (17)  for  C02(J"=1)  R  transition  perturbed  by  J-0  level 
of  N2. 


Figure  8.  Percent  difference  between  diffusion  operator  calculated  from 
equations  (16)  and  (17)  for  C02(J”=1)  R  transition  perturbed  by 
J=0  levei  of  N2. 
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moment  squared  from  the  line  intensity  for  most  of  the  transitions  on  the 
database.  In  addition  we  have  slightly  modified  the  parameter  on  the 
database  from  the  standard  definition  of  the  transition  moment  squared. 
This  is  to  facilitate  the  use  of  the  parameter  for  computation  of  Einstein 
coefficients.  Below  we  describe  the  relationships  between  the  line 
intensity,  Einstein  coefficients,  transition  moments  and  oscillator  strengths. 

We  start  by  defining  the  number  of  molecules  per  unit  volume  in 
quantum  states  /  (lower  state  labeled  by  vibration,  rotation  quantum 
numbers  n,J,M)  and  u  (upper  state  n’,J',M’) 


N/ 


=:NM 


•Ej/kT 


N« 


Ngwe 


-E./kT 


(18) 


where  N  is  the  total  number  of  molecules  per  unit  volume,  g/  and  gu  are 
the  degeneracies  of  states  /  and  u  respectively,  and  Q  is  the  total  internal 
partition  function.  We  consider  isotropic  space  and  non  degenerate  levels 
at  first  and  later  generalize  to  the  degenerate  case.  The  Einstein 
coefficients  are 

B/..>|4  Einstein  coefficient  for  induced  absorption 

Bu-*i  Einstein  coefficient  for  induced  emission 

Au~*i  Einstein  coefficient  for  spontaneous  emission. 

The  number  of  transitions  per  unit  time  per  unit  volume  from  state  / 
to  u  is 


N/B/_*upvi.  *  Nj4Bj4_»/pVfcj  (19) 

and  the  number  of  transitions  per  unit  time  per  unit  volume  from  state  u 
to  /  is 


(20) 


At  equilibrium  there  is  u  balance 

bl/B/_*upV(,  =  Ny  {By_,/Pv^j  +  .Ay-,/) 


(21 ) 
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Rearranging  equation  (21)  gives 


N/  _  ^u-*l 

Na  B/_>mPv,, 


which  is  well  known  from  Boltzmann  statistics 

Nj.,  Ng,c-W  Q-i 

Combining  equations  (22)  and  (23)  yields 

p  _ _ Aa-»//Ba-»/ _ 


Comparing  this  with  the  Planck  function, 

8^ - 1 - 

c3  -  1  . 


gives 


and 


8/  ®/-*U  “  8m  ®M-*i 

The  integrated  absorption  coefficient  can  be  written 

^/m  “  I  ^^L.vdv  =  (N/B/_»a  -  Nh  Ba_,/) 

J  Av/k 


(22) 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 
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and  using  Eq.  (27)  gives 


(29) 


or  from  Eq.  (26)  we  End 


8iiv]„  \  N/g«  / 


(30) 


At  thermal  equilibrium  we  have  from  equation  (6)  N«g//N/gu  =  exp(- 
hv/u/kT)  and  equations  (29)  and  (30)  become 


S/„  =  N/B/_Jl-eh'"^T)h^ 

c 


(31) 


Siu  =  ) 


87iv]„ 


The  quantity  S’/«  represents  the  integral  of  the  spectral  absorption 
coefficient  over  effective  width  Av/^  of  the  spectral  line  whose  center  is  at 
v/m.  It  has  dimensions  of  cm-lsec'l  and  is  referred  to  as  the  integrated 
absorption.  It  practice  one  uses  a  spectral  absorption  coefficient  in  enr 
>atm-*  given  by  S’/^/pc  (p  in  atms) 


S/„  =  N/B/_  jl  -  ehV'.AT' ) 

pc  c  pc 


S,„  =  Nl  A„^,  Ml  .  ek'-AT  ) 
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We  now  wish  to  relate  this  to  the  transition  moment  squared.  From 
Penner*5  we  find  the  transition  moment  squared  is  related  to  the 
probability  that  in  unit  time  a  system  will  undergo  a  transition  from  a 
state  n  to  a  state  m  under  the  influence  of  an  external  field  of  volume 
density  pvmn  where  Em>En  by 

t  3h^  (35) 

Relating  m->tt  and  n-»Z  the  probability  of  a  transition  is  also  given  by 


Cm  Cm 
t 


—  pv,, 


(36) 


From  equations  (26),  (27),  (35),  and  (36)  we  find 


3h2 

(37) 

3h^ 

(38) 

3hc3  8« 

(39) 

Finally,  from  equations  (33)  and  (37),  or  (34)  and  (39)  we  find 


3hc2  P 


(40) 


or  with  v/m/c  =  (o/m 


3hc  P 


(40a) 


Equations  (40)  relate  the  line  intensity,  S/u,  to  the  transition  moment 
squared,  IR/ul^,  for  the  non  degenerate  case. 

The  dimensionless  absorption  oscillator  strength  //-»u  is  defined 
through  the  relation 


where  e  and  m  denote,  respectively,  the  electronic  charge  (in  esu)  and  the 
mass  (in  gm/electron).  Hence  the  integrated  absorption  and  dimensionless 
/-number  are  related  through  the  expression 


mc2  P 


(42) 


The  emission  oscillator  strength  fu->l  is  related  to  //_»«  through  the 
expression 


fu-*l  — 

g« 


(43) 


From  equations  (37),  (40a),  and  (42)  we  can  relate  the  absorption  oscillator 
strength,  the  transition  moment  squared,  and  the  Einstein  coefficient  for 
induced  absorption 


3he2  ji2e2 


(44) 


Relationships  with  the  other  Einstein  coefficients  can  be  easily 
obtained  from  Eqs.  (38)  and  (39). 

States  are  usually  degenerate  and  we  must  add  this  complication  to 
the  formulas.  For  example,  evaluation  of  the  number  of  molecules  in  a 
particular  state  from  Boltzmann  statistics  in  the  most  general  case  is 

_  Nd.e-E-/kT  _  Ngigj(2J-H)(2S+l)(2-6oA)eW 

Q  Q  (43) 
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where  dj  is  the  degeneracy  of  state  i  given  by  the  coupling  of  the  nuclear 
spins  with  the  rotational  wave  functions,  gi,  the  state  independent  nuclear 
spin  factors,  gj=nk(2Iic+l),  the  M  degeneracy,  (2J+1),  the  spin  degeneracy, 
(2S+1),  and  lamda  doubling  (2-8oa)^  and  Q  is  the  total  internal  partition 
function. 

Let  us  start  by  considering  a  simple  case  with  M  degeneracy  only. 
The  formulas  derived  so  far  apply  to  the  non  degenerate  case  thus  we  can 
write  the  following  expressions  for  the  non-degenerate  sublevels: 

A  (nJMXn’JM*)  =  |R(nJM)(n'JM’f 

3hc3  .  (46) 

The  integrated  absorption  for  the  non-degenerate  case  is  then 

S(nJM)(n’JM')  =  ^  A(nJM)(n’JM*)  (l  -e**'""^'''*') 

P 

3hc  P  .  (47) 

However  in  the  absence  of  an  electric  or  magnetic  field,  the  levels  are 
degenerate  with  degeneracies  gnj  =  (2J+1)  for  the  lower  state  and  gnj*  = 
(2J’+1)  for  the  upper  state.  For  transitions  between  the  degenerate  energy 
levels  En’j*  and  Enj,  the  fundamental  expression  for  the  Einstein  coefficient 
for  spontaneous  emission  becomes 

A(nJ)(n'J)  =-^  2  A(nJM)(n'J-M') 

M.M'  ,  (48) 

where  the  summation  extends  over  all  allowed  values  of  M  and  M’  that  are 
consistent  with  the  specified  energy  transition.  Equation  (48)  follows  from 
the  fact  that,  in  the  absence  of  a  radiation  field,  the  Einstein  coefficient  is 
defined  in  such  a  way  that 


=  NnJM  A(nJM)(n'J'M')  =  A(nJM)(n'J'M') 
d  t  8nJ 


(49) 


whence 


.dN^-  = 

dt 


2  =  2  A(nJM)(n-JM-) 

M.M’  °  ^  M.M’ 


(30) 


Equations  (49)  and  (50)  show  that  Nn*j’  =  Nn’j’M'*(2J’+l),  i.e.  (2J+1)  times 
other  degeneracies.  Using  (46),  equation  (48)  can  be  rewritten  for  the 
degenerate  levels 


A(nJ)(n'J')= 


3 

(nJ)(n’J*) 


z 

M.M 


^^^|R(nJM)(n'J'M'f 
.  3hc3 


(51) 


Using  (47)  for  non-degenerate  energy  levels  along  with  (48)  we  find  that 

S(nJ)(n'J)  =  Z  S(nJM)(n'JM’) 

M.M-  (52) 


or 

S(nJ)(irj)  =  t*J(nJ)(n‘J)  (l  ’  J  |R(nJM)(n'JM')| 

3hc  P  gnJ^^. 

We  now  consider  the  general  situation  where  the  lower  and  upper  states 
may  have  several  types  of  degeneracy,  for  example  a  state  may  have  M 
degeneracy  from  the  various  projections  of  J  on  a  space  fixed  axis  and  a 
degeneracy  from  the  coupling  of  nuclear  spins  with  the  rotational 

wavefunctions.  The  quantum  numbers  labeling  the  states  can  be  divided 

into  two  sets;  one  set  for  the  non-degenerate  labels,  q,  and  the  other  set 
labeling  the  degeneracies.  C-  T'o  illustrate  this,  above  we  have  been 

considering  the  vibration-rotation  state  labeled  by  n.J.M.  In  this  labeling 

scheme  the  non-degenerate  set  of  quantum  numbers  ,q  ,  is  n  and  J,  and  the 
degenerate  set,  is  the  M’s.  The  energy  of  the  state  is  labeled  by  q  and 
the  number  of  molecules  in  each  state  is  given  by 
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The  Einstein  coefficient  for  spontaneous  emission  is  defined  from 


4.^- 


dt 


and 


dN 


^=2 


•dN^.4-  _  Nrj 


dt  rt.  dt  d 


2  A(ti4)(tt4-) 


(55) 

(56) 


where  the  summation  extends  over  all  allowed  values  of  that  are 
consistent  with  the  specified  energy  transition.  Again  we  note  that  Nt^  = 
dn*Nt^4  The  Einstein  coefficient  for  spontaneous  emission  becomes 


dn-  5.V 


(57) 


or 


Ari’-»n  =  2  |R(Ti^)(Ti'^’)f 

3  n  Qti* 


(58) 


We  also  have 


2  |R(Ti^)(ti’^')F 
3h^  dll'  4.V 


(59) 


and 

J-  2  |R(Tl4)(ll'd^ 

3h'  U-  .  (60) 


In  this  notation  the  line  intensity  is  related  to  the  transition  moment 
squared  by 


(61) 


Using  equations  (54)  and  (61)  we  find  the  following 


_L  J  =211£_E  _ Q(T) 

8k3  Nd,,  e-E-AT 


(62) 


where  S(n)(t^')  has  units  of  cm-Vatm.  Rewriting  (62)  in  HITRAN  units  of 
cm*V(molecules»cm'2)  p  and  N  are  absorbed  into  S  to  get 


JL  T  |R/„m/-.J2  =11ic.  _ Q(T^) 

8n3  (l- e^WkT ) 


(63) 


This  is  the  quantity  that  has  been  added  to  the  database  in  the  transition 
moment  squared  field  we  will  call  the  weighted  transition  moment 
squared.  This  is  more  useful  to  those  wanting  Einstein  coefficients  since 
this  is  just  3h2/2n  times  Bi_»u- 

The  program  TR-MOM.for  reads  each  line  in  HITRAN  format  and. 
depending  on  the  molecule,  decodes  the  necessary  information  (quantum 
numbers  etc.)  for  the  degeneracy  calculation,  calculates  the  total  internal 
partition  sum  (TIPS),  and  then  computes  the  weighted  transition  moment 
squared,  /?2.  The  line  is  then  written  to  a  file  with  this  additional 
information.  TR-MOM  has  been  run  on  all  the  lines  in  the  1991  HITRAN 
database.  For  molecules  for  which  we  do  not  have  the  TIPS  a  zero  is 
written  for  the  /?2 

There  are  several  advantages  to  the  definition  Eq.  (63).  First,  it 
removes  from  the  partition  functions  the  degeneracy  of  the  state  t).  This  is 
important  in  that  our  partition  functions  were  calculated  with  all  the 
degeneracy  factors  accounted  for.  Thus  our  Q(T)  could  be  used  to  calculate 
thermodynamic  quantities.  In  spectroscopic  applications,  the  state 
independent  degeneracy  factors  (e.g.  non-coupled  nuclear  spin)  are  often 
left  out.  Had  we  not  removed  the  factor  of  dj^  from  the  partition  sum,  the 
transition  moment  squared  on  the  database  would  have  carried  these 
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factors.  This  could  have  caused  many  problems  to  the  users.  Another 
reason  is  that  Rritj'.  as  we  define  it,  is  3h2/2jc  times  the  Einstein  coefficient 
for  induced  absorption. 

The  transition  moment  squared  can  be  expressed  in  terms  of  the 
matrix  elements  of  the  various  j  components  of  the  dipole  moment  by 


(64) 


where  n,  J,  and  M,  denote,  respectively,  the  vibrational,  rotational,  and 
magnetic  quantum  numbers;  |x’j(r,0,(j>)  is  the  effective  component  of  the 
electric  dipole  moment  vector  p(r,e,(j))  and  H'n.j.M.  stands  for  the  complete 
time-independent  wave  function  of  the  molecule.  We  assume  separability 
of  the  wave  function  in  the  form 


VnJM  ■  S„j(r)  gjM(0.<l>).  8jm(0*O)  =  ,  (65) 

the  j*!’  component  of  the  dipole  moment  p’j  may  be  written  as 

\ifrM)  =  lt(r)  fj(0,<|>)  ^  (55) 

with  p(r)  given  by 

i  .  (67) 


The  transition  moment  squared  can  now  be  written 

2  |R(nJM)(n'JM'f  =  |s„jJp(r)|Sn'j'f  -^2  2  |8j.M 

.M.  M'  M.  M'  j  .(68) 

For  molecules  with  vanishing  electronic  angular  momentum,  the  double 
sum  appearing  in  the  preceding  expression  is  the  reduced  dipole  moment 
matrix  element  for  the  rotational  transition  from  state  J  toJ’ 
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(69) 


^  IS  2  |gj.M  =(jll^llJ'f 

M.  M’  j 


These  matrix  elements  obey  the  sum  rule 

-1 

J*  .  (70) 

The  square  of  the  vibration-rotation  matrix  element  is  denoted  by 

=  |(s;.l»(r)|S„.jf  _  (7,) 


and  it  follows  that  with  the  definition  of  the  weighted  transition  moment 
squared  given  by  Eq.  (63)  and  if  the  vibration-rotation  matrix  element  has 
little  J  dependence  we  have 


2 

J’ 


— L, 

KM0.0F 


1 


(72) 


The  rotational  states.  J*  ,  to  sum  over  in  equation  (72)  are  determined  by 
the  selections  rules  for  the  particular  molecule  under  consideration. 

The  sum  rule  given  by  Eq.  (72)  allows  us  to  perform  checks  on  the 
line  intensities  and  weighted  transition  moment  squared  for  bands  on  the 
database.  Any  time  we  find  the  sum  not  being  close  to  one  we  should  be 
suspicious  of  the  data.  This  analysis  also  gives  the  progression  of  the 
vibration-rotation  matrix  element,  which  can  be  used  to  determine 
Hermann-Wallis  factors.  We  have  performed  this  check  on  some  of  the 
HITRAN91  data.  In  particular,  we  have  looked  at  CO2,,  CO.  H2O.  HF,  HCI,  and 
N2.  The  sum  on  all  connecting  states  depends  on  the  molecule,  for  example 
for  CO2  for  a  particular  J”  we  must  sum  on  the  P,  Q,  and  R  transitions  from 
that  level.  For  CO,  HF.and  HCI  only  P  and  R  branches  enter  into  the  sum;  for 
water  vapor  the  selection  rules  for  an  asymmetric  rotor  must  be  used. 
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From  the  following  relationship, 


X  lR(nJM)(n'JM')l^  =  X  |R(n’JM’)(nJM^ 

M.M’  M.  M*  ,  (73) 

one  can  show  that  dtj*RTiTi’=  dn’*RTi'T^.  Thus,  for  a  particular  vibrational 
band  if  one  has  the  value  of  the  reduced  dipole  matrix  element  from  state 
nJM  to  n’J'M’  and  the  value  for  n’J'M’  to  nJM  is  needed  one  can  uses  the 
formula 

(j'llulljf  =  ^(j|lnlljf 

dj'  .  (74) 

This  is  necessary  when  applying  the  sum  rule  to  water  vapor  for  example. 
It  is  also  useful  for  other  molecules  when  all  the  RT^n'’^  ^re  not  available. 
Table  8  gives  the  results  for  some  of  the  checks  done  using  Eq.  (72) 

From  these  studies,  we  have  found  a  typographical  error  for  the  R2 
line  of  the  fundamental  of  HF  (noted  above).  Comparing  the  line  intensities 
we  found  a  value  of  2.26  xlO**  from  the  1986  database  and  4.196  xlOi^^ 
from  the  1991  databases  After  discussing  this  with  R.  Tipping,  he  stated 
that  the  HF  data  was  created  by  hand  hence  typos  were  a  possibility. 
Tipping  thinks  that  the  4  should  be  a  2,  i.e.  the  1991  value  should  be  2.196 
xlQl^.  He  will  check  this  and  during  his  visit  to  the  Geophysics  Laboratory 
in  July  we  will  correct  the  data.  An  interesting  problem  was  found  in  the 
HCI  data  one  of  the  R^tj’  had  an  exponent  of  -61.  These  numbers  were 
created  by  the  FORTRAN  program  TR-MOM  with  no  hand  operations.  As  a 
check  the  HCI  data  was  extracted  and  TR-MOM  run  again  on  this  and  all  the 
data  was  correct.  Thus  we  are  not  certain  at  this  time  how  the  E-61 
occurred  on  the  database.  We  will  consider  writing  a  general  program  to 
check  the  Rtjh-  on  the  database. 


Table  8.  Sum  Rule  Checks  of  R 
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2  2  1 
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1  1  0 
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♦**  Error  in  line  intensity 


7.0  1991  HITRAN  WORKSHOP 

The  1991  HITRAN  workshop  was  held  at  the  Geophysics  Laboratory 
from  June  13'**  to  the  14‘l‘  There  were  many  interesting  talks  on  both  the 
status  of  the  1991  Edition  of  the  database  and  on  future  needs  and  data. 
Below  the  final  session  entitled  Summary,  Review,  and  Recommendations 
is  summarized,  molecule  by  molecule. 
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Few  comments  about  water  vapor,  V.  M.  Devi  mentioned 
having  some  shift  measurements. 

Both  V.  M.  Devi  and  A.  Barbe  said  that  they  had  pressure  shift 
measurements.  A  problem  with  the  CO2  halfwidths  on  the 
database  (Rosenmann  et.  al.*®**!)  being  too  large  was  brought 
up.  In  general  the  halfwidths  are  needed  to  better  than  5% 
uncertainty.  Rothman  mentioned  the  Riii^*’s  being  fully 
implemented.  J.  Johns  questioned  the  definition  we  u.sed  but 
when  he  found  it  is  related  to  Hertzberg’s'^*  definition  of  the 
Einstein  coefficient  for  induced  absorption  he  said  it  was  ok. 

Not  much  said. 

Toth  will  have  new  data  soon,  also  pressure  shifts  for  v3. 

Not  much  said,  L.  Brown  requested  a  copy  of  Tipping’s  final 
report  on  CO. 

Discussion  centered  around  the  unassigned  lines.  L.  Brown  set 
the  lower  state  energy  to  300.555  whereas  HITRAN  has  used  a 
negative  E”  in  the  past.  It  was  suggested  to  change  Linda’s 
values  to  negative  numbers.  M.  A.  H.  Smith  has  some  CH3D 
data. 

Need  a  complete  revamping. 

No  comments. 

Update  needed,  A.  Pine  has  V1+V3. 

Use  data  of  Perrin  et.  al.. 

Update  with  data  of  S.  Urban.  GEISA  has  data  from  Meudon 
with  many  unassigned  lines.  There  are  notation  problems  -  not 
consistent. 

No  comments. 

Goorvitch,  Chackerian  data.  J.  Selby  questioned  A.  Goldman  on 
the  intensities  being  off  by  a  factor  of  2  to  6. 

These  are  fine. 

Randy  May  to  measure  Hermann-Wallis  factors.  K.  Chance  said 
it  would  be  nice  to  have  cross  ladder  transitions.  A,  Goldman 
mentioned  he  may  have  them. 

Possibly  an  update. 

N.  Husson  has  data. 
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H2CD  No  comments. 

H2O2  Questions  on  the  Hillman  data,  J.  Johns  has  made 
measurements,  A.  Perrin  also  has  data. 

C2H2  No  Comments. 

N2  No  Comments. 

PH3  New  data  available  from  Tarrago. 

SF6  No  comments. 

HCOOH  No  Comments. 

H2S  No  Comments. 


J.  Selby  commented  on  TIPS.  He  would  like  to  see  the  high 
temperature  partition  sums  improved,  especially  from  2000-3000K. 
C.  Benner  commented  on  reference  and  error  codes.  They  are 
extremely  useful  and  should  always  be  implemented. 

S.  Clough  wanted  to  know  when  the  next  edition  would  be 
distributed.  This  created  a  discussion  which  indicated  that  the 
scientific  community  wants  yearly  updates  if  possible. 

Other  comments  dealt  with  the  addition  of  coupling  coefficients  to 
the  database  in  the  future.  No  consensus  was  arrived  at.  Here  I  would  like 
to  add  some  comments.  At  this  phase  of  evolution  of  the  database  1  think 
it  might  be  useful  to  add  fields  for  the  temperature  exponent  for  self¬ 
broadening,  as  well  as  error  and  reference  codes  for  self-broadening  and 
shift  parameters.  However  I  do  .,ot  believe  that  adding  coupling 
coefficients  to  the  database  would  be  wise.  All  the  spectroscopic  data  that 
are  presently  on  the  database  are  fixed  by  nature.  By  this  1  mean  that  if 
se\era!  laboratories  make  measurements  or  calculations  are  made  of  the 
same  quantity  they  should  agree,  i.e.  there  is  one  "God”  given  value,  if  you 
like.  With  the  coupling  coefficients  there  is  a  danger  in  that  they  are 
dependent  on  the  model  chosen.  Even  worse,  within  the  model  the 
coupling  coefficients  are  determined  via  sum  rules  over  a  number  of  lines 
and  when  applied  all  the  lines  must  be  used  or  serious  errors  do  occur  as 
shown  by  L.  Strow.  It  would  be  dangerous  to  mix  these  two  types  of  data. 
1  would  suggest  the  inclusion  of  coupling  coefficients  and  any  other  of  this 
type  of  data  as  auxiliary  files  added  to  the  third  file  of  the  database.  The 
problem  with  adding  this  type  of  data  to  the  main  database  is  that  users 
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take  this  data  as  “God”  given  and  believe  fully  in  results  generated  from 
them. 


8.0  PAPERS  AND  PRESENTATIONS 

8.1  The  43‘h  Symposium  on  Molecular  Spectroscopy,  The  Ohio  State 
University.  June  13*17,  1988. 

Paper  Tell. ‘Simplified  Exprjssion  for  the  Total  Internal  Partition 
Sum  as  a  Function  of  Temperature”  with  Robert  Hawkins  and  Laurence  S. 
Rothman. 

8.2  1989  HITRAN  Workshop 

On  June  8-9,  1989  a  workshop  on  the  HITRAN  database  was  held  at 
the  Geophysical  Laboratory  at  Hanscom  Field.  I  gave  two  talks  at  the 
workshop,  the  first  reviewed  the  collision  broadened  halfwidths  and  their 
temperature  dependences  and  the  second  discussed  the  program  SELECT 
along  with  some  aspects  of  the  database.  The  first  talk  emphasized  the 
updates  expected  for  the  air*,  and  self-broadened  halfwidths,  the  line 
shifts  and  the  temperature  exponent  of  the  air-broadened  halfwidth.  The 
talk  was  to  allow  the  audience  to  have  an  opportunity  to  request  certain 
data  to  be  put  on  for  these  fields,  to  bring  to  our  attention  data  that  was 
overlooked  in  our  search,  and  to  bring  to  our  attention  unpublished  data. 
John  Ballard  of  Rutherford  Appleton  Laboratory  told  us  about  some  water 
vapor  data  of  his  student  John  Remedies  which  is  unpublished  and  about 
some  recently  published  NO  data.  John  will  send  the  data  to  us  in  July. 
Kelly  Chance  of  Harvard-Smithsonian  Center  for  Astrophysics  has  sent  me 
his  results  for  OH  and  will  keep  me  informed  of  some  ongoing 
determinations  of  broadening  parameters  for  this  molecule. 

The  second  talk  was  on  SELECT  and  presented  the  new  features  and 
options  of  the  program.  These  include  options  for  increased  speed  for 
searches,  version  labels  that  will  link  SELECT  and  HITRAN,  output  at 
temperatures  other  than  the  reference,  296  K,  and  more.  The  current 
error  codes  were  shown  and  were  discussed  with  the  audience.  This  led  to 
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the  codes  for  the  line  intensity  and  the  halfwidths  being  changed.  The 
codes  now  follow  the  wavenumber  form  in  that  as  the  code  number  gets 
larger  the  accuracy  of  the  data  'ncreases.  It  was  also  decided  to  keep  the 
error  codes  for  S  and  y  the  same.  The  implementation  of  the  reference 
code  and  methods  for  its  use  was  also  presented  and  input  from  the 
community  requested. 

8.3  Journal  of  Molecular  Spectroscopy 

R.  R.  Gamache,  R.  L.  Hawkins,  and  L.  S.  Rothman,  “Total  Internal 
Partition  Sums  in  the  Temperature  Range  70-3000  K:  Atmospheric  Linear 
Molecules,”  J.  Mol.  Spectrosc.  142,  205-219  (1990). 
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